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ABSTRACT. Aperiodic tilings of Euclidean space can profitably be studied from 
the point of view of dynamical systems theory. This study takes place via a 
kind of dynamical system called a tiling dynamical system. 


1. Introduction 


In this chapter we study tilings of Euclidean space from the point of view of dy- 
namical systems theory, and in particular, symbolic dynamics. Our goal is to show 
that these two subjects share many common themes and that they can make useful 
contributions to each other. The tilings we study are tilings of Rt by translations 
of a finite number of basic tile types called “prototiles”. A good general reference 
on tilings is [@S87]. The link between tilings and dynamics will be established 
using a kind of dynamical system called a tiling dynamical system, first described 
by Dan Rudolph [Rud88], [Rud89]. The parts of the theory we concentrate on 
here are the parts most closely related to symbolic dynamics. Interestingly, these 
also tend to be the parts related to the theory of quasicrystals. 

A quasicrystal is a solid which, like a crystal, has a regular enough atomic struc- 
ture to produce sharp spots in its X-ray diffraction patterns, but unlike a crystal, 
has an aperiodic atomic structure. Because of this aperiodicity, quasicrystals can 
have “symmetries” forbidden to ordinary crystals, and these can be observed in 
their X-ray diffraction patterns. The first quasicrystals were discovered in 1984 by 
physicists at NIST (see [SBGC84]) who observed a diffraction pattern with 5-fold 
rotational symmetry. For a good mathematical introduction to quasicrystals see 
[Sen95]. 

The theory of quasicrystals is tied up with some earlier work on tiling problems 
in mathematical logic ([Wan61], [Brg66]). Central to this circle of ideas is the 
concept of an aperiodic set of prototiles. One of the most interesting aperiodic sets, 
which anticipated the discovery of quasicrystals, is the set of Penrose tiles, discov- 
ered in the early 1970s by Roger Penrose [Pen74]. Penrose tilings play a central 
role in the theory of tiling dynamical systems because they lie at the crossroad of 
the three main methods for constructing examples: local matching rules, tiling sub- 
stitutions, and the projection method. As we will see, the tiling spaces constructed 
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by these methods are analogous to three well known types of symbolic dynamical 
systems: finite type shifts, substitution systems and Sturmian systems. 

The connection between tilings and symbolic dynamics goes beyond the analo- 
gies discussed above. Since tilings are (typically) multi-dimensional, tiling dynamics 
is part of the theory of multi-dimensional dynamical systems. We will show below 
that one can embed the entire theory of Z? symbolic dynamics (the subject of Doug 
Lind’s chapter in this volume) into the theory of tiling dynamical systems. It turns 
out that much of the complication inherent in multi-dimensional symbolic dynamics 
(what Lind calls “the swamp of undecidability” ) is closely related to the existence 
of aperiodic prototile sets. 

Finally, one can view tiling dynamical systems as a new type of symbolic dy- 
namical system. Since tilings are geometric objects, the groups that act naturally 
on them are continuous rather than discrete (i.e., R? versus Zt). Because of this, 
one needs to define a new kind of compact metric space to replace the shift spaces 
studied in classical symbolic dynamics. We call this space a tiling space. Even in 
the one dimensional case (i.e., for flows) tiling spaces provide a new point of view. 

In the first part of this chapter we carefully set up the basic theory of tiling 
dynamical systems and give complete proofs of the main results. In later sections, 
we switch to survey mode, giving references to access the relevant literature. Of 
course there are many topics we can not cover in such a short chapter. 

These notes are based on a AMS Short Course presented by the author at the 
2002 Joint Mathematics Meeting in San Diego, California. The author wishes to 
thank Tsuda College in Tokyo, Japan and the University of Utrecht, The Nether- 
lands, where earlier versions of this course were presented. My thanks to the Natalie 
Priebe Frank and Cliff Hansen for carefully reading the manuscript and making 
several helpful suggestions. My thanks also to the referee who suggested several 
substantial improvements. 


2. Basic definitions in tiling theory 


2.1. Tiles and tilings. A set D C Rf, d > 1, is called a tile if it is com- 
pact and equal to the closure of its interior. We will always assume that tiles are 
homeomorphic to topological balls, although in some situations it is useful to allow 
disconnected tiles. Tiles in R are closed intervals. Tiles in R? are often polygons, 
but fractal tiles also occur frequently in examples. 

A tiling! x of R? is a collection of tiles that pack Rê (any two tiles have pairwise 
disjoint interiors) and that cover R? (their union is R¢). Two tiles Dı, D2 are 
equivalent, denoted Dı ~ Dz, if one is a translation of the other. Equivalence class 
representatives are called prototiles. 


DEFINITION 2.1. Let T be a finite set of inequivalent prototiles in R°. Let Xy 
be the set of all tilings of R? by translations of the prototiles in 7. We refer to Xz 
as a full tiling space. 


Broadly speaking, geometry is concerned with properties of objects that are 
invariant under congruence. Similarly, dynamics is generally concerned with group 
actions. In this chapter, we will be interested in how groups of rigid motions act on 


lWe use the lower case notation x for a tiling because we want to think of x as a point ina 
tiling space X. 
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sets of tilings. Because of this, we will distinguish between congruent tilings in Xy 
that sit differently in Rt. Of central interest will be the action of R? by translation. 


DEFINITION 2.2. For t € R? and x € Xz let Ttr € Xz be the tiling of Rt in 
which each tile D € x has been shifted by the vector —t, that is Ttx = {D — t : 
D € x}. We denote this translation action of R? on Xz by T. 


The primary reason for studying T is that it is related to the long range order 
properties of the tilings in Xz. While such properties are geometric in nature, we 
will gain access to them through dynamical systems theory. 


2.2. Local finiteness. Let 7 be a set of prototiles. A T-patch y is a finite 
subset y C x of a tiling x € Xz such that the union of tiles in y is connected. This 
union is called the support of y and written supp(y). The notion of equivalence 
extends to patches, and a set of equivalence class representatives of patches is 
denoted by 7*. The subset of patches of n tiles, called the n-patches, is denoted 
by T™) C T*. 

We will impose one additional condition, called the local finiteness condition, 
on all tiling spaces X7. 


DEFINITION 2.3. A tiling space Xy has finite local complexity if T® is finite. 


Equivalently, 7”) is finite for each n. Sometimes the geometry of the tiles 
themselves will impose the local finiteness condition, but we usually need to add 
it as an extra assumption. From now on, whenever we write 7, 7* or Xz, it will 
always implicitly include a choice of a finite T°). When working with polygonal 
prototiles in R?, a common way to achieve local finiteness to assume that all tiles 
meet edge-to-edge. 


EXAMPLE 2.4. Consider the set S consisting of a single 1 x 1 square prototile. 
Without any local finiteness condition, fault lines exist in the tilings x € Xs with 
a continuum of possible displacements. Imposition of the edge-to-edge condition 
means that every x € Xs is a translation of a single periodic tiling. See Figure 1. 


(a) (b) -1 (c) 


FIGURE 1. (a) Part of an edge-to-edge square tiling. (b) A square 
tiling with a fault having displacement t. (c) Local finiteness can 
always be forced geometrically by cutting “keys” on the edges of 
tiles. 


EXAMPLE 2.5. We get more interesting square tiling examples by taking Sn, 
n > 1, to be the set of 1 x 1 square prototiles marked with “colors” 1,2,...,n. To 
do this we also need to modify our notion of equivalence so that differently colored 
squares are not considered to be equivalent. 
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Now consider the subset Xo C Xg, consisting of all tilings whose vertices lie 
on the lattice Z” C R? and let To be the restriction of the R? shift action T to the 
subgroup Z4. It is clear that Xo is Tp-invariant. 


We will see later how this example links tilings to discrete symbolic dynamics. 


EXAMPLE 2.6. Fix n > 4. Let s = n for n odd, and s = 2n for n even. For 
0<k <n let vk = (cos(27k/s),sin(27k/s)) € R?, i.e., vk is a sth root of unity, 
viewed as a vector in R?. Let R,, denote the set of all > rhombi with translations 
of the vectors v as sides. Define Xr, to be the corresponding edge-to-edge tiling 
space. Two examples of x € Xr, (one with markings) are shown in Figures 3 and 
10. 


2.3. The tiling topology. As we now show, finite local complexity tiling 
spaces have particularly nice topological properties. The tiling topology is based 
on a simple idea: two tilings are close if after a small translation they agree on a 
large ball around the origin (see [Rad99], [Rob96b], [Rud88], [S0197]). However, 
the details turn out to be a little subtle. 

Given K C R? compact and x € Xz, let x|[K]] denote set of all sub-patches 
xz’ C x such that K C supp(a’). The smallest such patch is denoted z|K]. For 
r > 0 let B, = {t € R° : ||t|| < r}, where ||- || denotes the Euclidean norm on Rt. 


LEMMA 2.7. For x,y € Xr define 


(2.1) d(x,y) = inf ({V2/2}U {0 < r < V2/2:42’ € 2[[B,/,]], 
y' € y[[Bi/,]], with T*x’ = y' for some ||t|| < r}). 
Then d defines a metric on Xr. 


We call d the tiling metric. 


PROOF. We prove only the triangle inequality. Let 0 < d(x,y) =a’ < d(y, z) = 
b with a’ +b! < V2/2. Let 0 < e < V2/2— (a! +0’) and put a = a’ + €/2 and 
b = b' + 6/2. Then there are 2’ € 2[[Byyal], y € yl[[Biyall, y” € yl[[Biys]] and 
2" € z[[B1,,]], and also t,s € R? with ||t|| < a and ||s|| < b, such that Ta’ = y' 
and T782z" = y”. 

Let yo = y' N y”, zo =T*tyo C x’ and zo = TSyo C z’. Then 


(2.2) T t+) z) = zo where ||t +s|| < a +b. 


Letting c = a + b, then since 0 < a < b < v2/2, 


E es 
>e Ge 
and it follows that By;, C (By, +t). Now y’, y” € y[[B1/b]] so xo € 2[[Bi, +t] € 


z[[Bi;d]]. 
Combining this with (2.2), we have d(x,z) < a+b = d(x, y)+d(y,z) +€, where 
€ > 0 is arbitrarily small. The triangle inequality follows. 


LEMMA 2.8. The tiling metric d is complete. 
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PROOF. Consider a Cauchy sequence £n of tilings. Assume d(£n+1, £n) > 0 
and let sn = d(£n41, £n) +27". By passing to a subsequence, we may assume 
Sn is decreasing and X] Sn < oo. it follows from (2.1) that for each n there 
exists tn € R? with ||t,|| < sn and z’, € £n[[B1/s„]] such that Tt» z!, C 2}... Put 

co 
rn = > py tk- Then 


r J Tnit / Patil 
Pet, = TT ae, CT ey 


This implies that T™2/, is an increasing sequence of patches, so we can define a 
tiling x = UnT™” x}. Finally, d(x, £n) < max(||rn||,5n) > 0. 


THEOREM 2.9. (Rudolph [Rud89]) Suppose Xr is a finite local complexity 
tiling space. Then Xz is compact in the tiling metric d. Moreover, the action T of 
Rê by translation is on Xz is continuous. 


EXERCISE 1. The proof of compactness amounts to the observation that the 
local finiteness condition is equivalent to Xy being totally bounded (see [Mun75)). 
Fill in the details of this proof and also prove the continuity of T. 


EXERCISE 2. Prove the following converse to Theorem 2.9: a translation in- 
variant set X of tilings which is compact in the tiling metric (2.1) must have finite 
local complexity. 


3. Tiling dynamical systems 


3.1. Tiling spaces as symbolic dynamical systems. Throughout this 
chapter, a dynamical system will be a pair (X,T) where X is a compact metric 
space (the phase space) and T is a continuous action of a group, usually (but not 
always) R. The study of the topological properties of dynamical systems is called 
topological dynamics. The study of the “statistical properties” of dynamical systems 
is called ergodic theory. An excellent introduction to both topological dynamics and 
ergodic theory is Walters? [Wal82]. 

Symbolic dynamics studies a special kind of dynamical system called a symbolic 
dynamical system. The classical set-up is 1-dimensional, but we describe here the 
general d-dimensional case (see also the chapter by D. Lind). For the group we take 
Zt, and we let X„ = {1,.. ny}, n > 1, with the product topology. Letting T 
be the shift action of Z? on Xn, we obtain a dynamical system (Xn, T) called the 
d-dimensional full shift on n symbols. In some ways this example itself is too simple 
to be interesting, but it has very complicated subsets. A Z4-symbolic dynamical 
system is defined to be a pair (X, T) where X is a closed T-invariant subset X C Xn 
called a shift space? (see [LM95] and [Que87}). 


DEFINITION 3.1. Let Xy be a full d-dimensional tiling space and let T denote 
the translation action of R?. A tiling space X is a closed T-invariant subset X C X7. 
We call the pair (X, T) a tiling dynamical system. 


Now we can precisely state our way of thinking of tiling dynamical systems as 
a new type of symbolic dynamical system. We think of the prototiles D € T as the 
symbols. The full tiling space Xy corresponds to the full shift, and more general 
tiling spaces correspond to more general shift spaces. Like the product topology, 


2Even though this book concentrates almost exclusively on Z actions, the theory goes through 
with very little effort to actions of Z4 and R4. 
3A shift space is also sometimes called a subshift. 
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the tiling topology is compact and metric, and in both cases closeness corresponds 
to a good match near the origin. However, in the the case of tilings there is the 
possibility of a small translation, and since we want this to be continuous, the 
topology needs to be defined accordingly. 


REMARK 3.2. The theory of tiling dynamical systems contains the theory of Z4 
symbolic dynamics. One can show that the space Xo constructed in Example 2.5 
is homeomorphic to the symbolic full shift, and that To implements the shift action 
on Xo. 


3.2. Finite type. Let Xz be a full tiling space and let F C T*. Let X\¢ C 
Xy be the set of all tilings x € Xz such that no patch y in x is equivalent to any 
patch in F. We call such a set F a set of forbidden patches. 

One can show that for any F C T*, the set X\¢ is a tiling space (i.e., it is 
closed and T-invariant). Moreover, it is clear that every tiling space X C Xz is 
defined by a set F of forbidden patches. However, the set F is not unique! 


EXERCISE 3. Prove the three statements in the previous paragraph. Hint: See 
[LM95]. 


In symbolic dynamics, the most important kind of shift space is a finite type 
shift. The following definition introduces the corresponding idea in tiling theory. 


DEFINITION 3.3. A tiling space X C Xz is called a finite type tiling space if 
there exists a finite F C T* so that X = X\F 


The most common case is F C 7®. This is called a local matching rule. It 
is convenient to formulate this case in terms of the allowed 2-patches rather than 
forbidden ones. To accomplish this, we put Q = T, let Q?) = T)\F, and write 
Xo for Xy. Note that imposing a local matching rule really just amounts to 
strengthening the local finiteness condition. Thus a full tiling space is a kind of 
finite type tiling space. 


EXAMPLE 3.4. (The Penrose tiles) Consider the marked version P (shown in 
Figure 2) of the prototiles Rs. The set P) (which defines the matching rules) 


Lf Æ 


FIGURE 2. The Penrose tiles. The protoset P consists of the two 
marked tiles shown, and all rotations so that edges have have angles 
27n/10. In particular, card(P) = 20 


imposes the requirement that the markings on any pair of adjacent tiles must match. 
(As we will see below, markings are often used for this purpose). We call P, together 
with the matching rules, the Penrose tiles. Tilings x € Xp are called Penrose 
tilings. Part of a Penrose tiling is shown in Figure 3. 
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FIGURE 3. A part of a Penrose tiling. 


3.3. The Tiling Problem. Suppose we are given a set 7 of prototiles and a 
set F C T* of forbidden patches. Consider the following problem: 


TILING PROBLEM. Is X\¢ # 0? 


We begin with a positive result, a version of which appeared in [Wan61] (see 
[GS87] for a proof). 


EXTENSION THEOREM. Let T be a collection of prototiles with a local finiteness 
condition T°?) and let F C T* be a set of forbidden patches. Define T+ C T* to 
be the set of patches that do not contain any forbidden sub-patches. Then X\ r #0 
if and only for each r > 0 there is a patch y € T* with B, +t C supply) for some 
teR?. 


The trouble with the Extension Theorem is that it is not constructive. To 
conclude that X\ z is nonempty one needs to see infinitely many patches in 7 +, 

This difficulty can be appreciated if one tries to tile the plane manually with 
Penrose tiles. There are a lot of “dead ends”: patches in y € P* that do not belong 
to P*. How can we know that P+ doesn’t have some largest patch y? Later, we 
will give a proof that Xp Æ 0, but that proof will require a new idea. 

The question of whether the Tiling Problem is decidable was raised by Wang 
[Wan61] for the case of marked square tiles Sn, together with a local matching 
rule. These are now known as Wang tiles. 
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When d = 1 there is an easy algorithm to answer the Tiling Problem: First we 
draw a graph G with vertex set 7 and directed edges 7), and let A be the m x m 
adjacency matrix for G, where m = card(T). The entries of A* give the number 
of paths of length k in G. If A™+! Æ 0 (i.e., not identically zero) then A* Æ 0 for 
any k > m, and we conclude Xo £ 0. 


DEFINITION 3.5. A tiling x of R? is called a periodic tiling if its translation 
group T, = {t € R? : Ttx = x} is a lattice: that is a subgroup of R? with d linear 
independent generators. A tiling x is called aperiodic if Ty = {0}. 


In the case d = 1 one can easily show that if Xy 4 @ then there is a periodic 
tiling x € Xz. Wang conjectured [Wan61] that the same holds for d > 1. 


WANG’S CONJECTURE. 


(1) There is an algorithm to decide the tiling problem. 
(2) Whenever Xr #9, there exists a periodic tiling x € Xr. 


Wang proved that (2) implies (1). The argument goes as follows: 

First suppose Xz Æ Ø. Then there exists a periodic x € Xz. For each n € N, 
list all tiling patches y € T+ with Ban C supp(y) and B,+41 Z supp(y). We will 
eventually see a complete period of x. In this case the algorithm will stop and 
answer “yes”. 

Now suppose Xy = Ø. Then it follows from the Extension Theorem that we 
will eventually find n € N so that no y € T* has support containing Bn. In this 
case the algorithm will stop and answer “no”. 

But Wang’s conjecture turns out to be false! For d > 2 Berger [Brg66] showed 
that the tiling problem is, in fact, undecidable. His solution included the construc- 
tion of an example of a 2-dimensional finite type tiling space Xg containing no 
periodic tilings. 


DEFINITION 3.6. 
(1) A nonempty tiling space X is called an aperiodic tiling space if it contains 
no periodic tilings (i.e., every x € X is an aperiodic tiling). 
(2) A prototile set Q with local matching rule Q®) is called an aperiodic 
prototile set if Xo is an aperiodic tiling space.* 


Berger’s original aperiodic prototile set Q satisfies card(Q) > 50,000 (see 
[GS87]). Later, Raphael Robinson [rRob71] found a simple example with card(Q) = 
32 (a picture of these tiles appears in Lind’s chapter of this volume). 


3.4. Counting prototiles and the “einstein” problems. The problem of 
finding small aperiodic sets of prototiles has been a popular one (see [GS87] for the 
history up to 1987). The exact formulation depends on how one counts prototiles. 
With our notion of equivalence (translation but not more general congruences), 
the Penrose tiles P consist of 20 prototiles. Counting this way, the current best 
example in R? is a set K of Wang tiles due to Kari and Culik ([Kar96], [Cul96}) 
with card(X) = 13. This example is particularly interesting because it is not related 
to any other known examples. 


4In some literature, the term “aperiodic tiling” is reserved for tilings x E€ Xg, where Q is an 
aperiodic prototile set. 
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It is perhaps more natural to allow congruence classes of prototiles to count 
only once (i.e., to allow rotations of the prototiles). With this system of counting 
there are just 2 Penrose tiles, and for R? this is the best result so far. 

The question of whether there exists an aperiodic prototile set consisting of a 
single tile has been named the “einstein problem” by Ludwig Danzer. For d = 3 
an example of an einstein was discovered by Schmitt and Conway (see [Sen95]). 
However, it tiles in a way that is, in a certain sense, very weakly aperiodic. 

There are some partial results on the 2-dimensional einstein problem as well. 
For any prototile set 7 = {D}, where D is a topological disk, whenever Xy 4 0 
there exists a periodic tiling x € Xy ([Ken92], [Ken93], [GBN89]). In other 
words, there is no einstein up to translation. 


4. Substitution tiling spaces 


4.1. Perfect decompositions. Let L € Gi(d,R) be an expansive linear trans- 
formation of R?. Expansive means that every eigenvalue of L lies outside the unit 
circle. The case L = AM, where M is an isometry and A > 1 is called a similarity. 

A perfect decomposition (or just a decomposition) is a mapping C : T — L~!T* 
that (up to equivalence) satisfies the perfect overlap condition 


(4.1) supp(C(D)) = supp(D), 

In a slight abuse of language, a decomposition is called self-similar if L is a 
similarity; in the general case it is called a self-affine decomposition. 

The mapping S = LC is called a perfect self-similar or self-affine tiling substi- 
tution” on T. 

In the case d = 1 there is no difference between a self-similar and a self-affine 
substitution, and any tiling substitution can be written S = AC, where A is a 
positive real number. In the self-similar case when d = 2, we can identify R? with 
C. By replacing S with S? we can assume L orientation preserving. Then we can 
regard L as multiplication by à € C, with |A| > 1, so that S = AC. 


EXAMPLE 4.1. (Polyomino decompositions) The chair decomposition Ce on 
the set C of four chair prototiles is obtained by taking the decomposition pictured 
in Figure 4(a) and its four rotations. The table decomposition C4, Figure 4(b), is 


fsa 


(a) (b) (c) 


FIGURE 4. (a) The chair, (b) the table and (c) the 3-dimensional table. 


defined on the protoset Də of two “dimers” in the plane. The 3-dimensional table, 
Figure 4(c), is defined on the set D3 of 6 dimers in R3. The asymmetry of this 


5This is also sometimes called an inflation mapping. 
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example makes it necessary to keep track of the prototiles’ orientations. Many 
other polyomino examples are easily devised. 
In all cases shown in Figure 4, L is a similarity with M = Id and \ = 2. A 


non-self-similar polyomino decomposition with L = (; J is shown in Figure 5. 


A EE ie 


FIGURE 5. The folding table: a non-self-similar version of the ta- 
ble. Here we show the decomposition of [D2 into Do. 


EXAMPLE 4.2 (Raphael Robinson’s triangular Penrose tilings). This is a non- 
polyomino self-similar decomposition C; with \ = (1/2)(1+ v5) on a set Pı of 40 
marked triangular prototiles. Two of the prototiles are shown in Figure 6(a). The 
decomposition Cı is shown in Figure 6(b). Let Xp, denote the finite type tiling 


AÀA 


FIGURE 6. The triangular Penrose tilings: The protoset Pı con- 
sists of a finite set of rotations of the two tiles shown that is closed 
under decomposition. 


space corresponding to the usual matching rule that the arrows on adjacent edges 
must match. We will show below how to use the tiling substitution S4 = AC}. to 
prove that Xp, # 0. 


EXAMPLE 4.3 (The pinwheel tilings). The self-similar decomposition shown in 
Figure 7 has been studied extensively by Radin (see [Rad94]). Up to rotation and 
reflection it has a single prototile. ie the decomposition, one copy of the prototile 
is rotated by an angle 0 = arctan D so that £ is irrational. Such a rotation 
has infinite order and hence there is no finite prototile set invariant under this 
decomposition. We can get around this difficulty by modifying the definition of 


equivalence and the tiling metric d to allow rotations as well as translations. 


The next two examples don’t quite satisfy (4.1), but can nevertheless easily be 
accommodated. We refer to them as imperfect decompositions and the correspond- 
ing substitutions as imperfect tiling substitutions (see Definition 5.25). 


EXAMPLE 4.4. (The Penrose decomposition). This decomposition applies to 
the marked prototile set P of Example 3.4. Since C(D) D D, this is an imperfect 
decomposition. 
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FIGURE 7. The pinwheel decomposition: One prototile is irra- 
tionally rotated relative to the others. 


a oe 


FIGURE 8. The imperfect decomposition of the rhombic Penrose tilings. 


EXAMPLE 4.5. (The binary tiling decomposition). The decomposition shown 
in Figure 9 is from [GL92]. A patch of a tiling x in the corresponding substitution 


LY LF 


FIGURE 9. The imperfect “binary” tiling decomposition C% is self- 
similar. In the corresponding tiling substitution Sẹ = AMC, the 
similarity L = AM includes a rotation M by 27/40. Note, however, 
that C?(T) C L~°T, so there is a version of S? with no rotation. 


tiling space is shown in Figure 10. A perfect version of this decomposition is shown 
in Figure 13. 
4.2. Properties of tiling substitutions. 


DEFINITION 4.6. Let S = LC be a tiling substitution on 7 where Xy has finite 
local complexity. We say S satisfies the 2-patch closure property if S(T) C T*. 


From now on, unless we say otherwise, we will assume every tiling substitution 
S satisfies the 2-patch closure property. The reason® for this assumption is that a 


6The need for such an assumption was pointed out to the author by Natalie Priebe Frank. 
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FIGURE 10. A patch of binary tiling. 


tiling substitution S : T — T* (satisfying 2-patch closure) can be extended to a 
mapping S : Xr — Xz. This mapping performs the decomposition C, viewed as a 
mapping C : Xr — X,-17, and then applies the linear expansion L to the entire 
tiling. Note that C satisfies translation invariance 


(4.2) TC = CT* for all t € R°. 


EXERCISE 4. Show that both C and S are continuous in the tiling topology 
and that S satisfies Tt S = ST*. 


We call a tiling substitution S = LC invertible if it is 1:1. In this case § has 
a continuous inverse ST! = L~'C~! on S(Xr). Then C7! is a continuous and 
translation invariant mapping, which is called a composition. 

Let us denote T = {Dj,..., Dn}. We define the structure matrix A of the tiling 
substitution S = LC to be n x n matrix with entries A; į equal to the number of 
prototiles equivalent to D; that occur in $(D,). A tiling substitution is called 
primitive if A! > 0 for some k > 0. 


4.3. Substitution tiling spaces. Let S be a tiling substitution on T. Take 
D € T and define a sequence of patches x, inductively: xı = {D}, and a, = 
S(xp-1), k > 1. Since we assume S' satisfies the 2-patch closure property, it follows 
that x, € T* for all k. 

The patches x, are used to define a tiling space as follows. First we define 
a set Fs of forbidden patches by stipulating that y € T* is forbidden if it is not 
a sub-patch of 2, for any k. The tiling space Xs = Xs, C Xv is called the 
substitution tiling space corresponding to S. The dynamical system (Xs,T) is 
called a substitution tiling dynamical system. 
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FIGURE 11. The patches zk, k = 1,2,3, for the chair tiling sub- 
stitution (Figure 4(a)), and a patch of chair tiling. 


LEMMA 4.7. Let Xs C Xy be the substitution tiling space corresponding to a 
primitive tiling substitution S on T. Then Xs #0 and Xs is independent of the 
initial tile D € T. Moreover S(Xg) C Xs. 


PROOF. First we observe that 2, € T+ for all k, since by the definition of 
Fs, no patch £k contains any forbidden sub-patches. We also note that since L is 
expansive, diam(supp(x,)) — oo. Thus given r > 0, there exists t € R? so that 
B, = t C supp(z,) for k sufficiently large. It follows from the Extension Theorem 
that Xs = Xr, #4. The independence of Xs from the choice of D follows from 
primitivity. The S-invariance is clear. 


COROLLARY 4.8. Let P be the set of marked rhombic Penrose tiles shown in 
Figure 2, and let Pı be the marked triangular Penrose tiles shown in Figure 6. Then 
the corresponding finite type tiling spaces Xp and Xp, are nonempty, i.e., Penrose 
tilings exist. 


PROOF. For the two tiling substitutions, S from Figure 8 and the substitution 
Sı from Figure 6, the corresponding substitution tiling spaces satisfy Xp D Xs 40 
and Xp, 2 Xs, F 0. 


This is essentially Penrose’s argument, although he did not couch it in dynami- 
cal terms. We conclude with a result that characterizes exactly which tilings belong 
to a substitution tiling space. 


PROPOSITION 4.9. Let Xs C Xz be a substitution tiling space for a primitive 
tiling substitution S. Let x € Xr. Then x € Xs if and only if there is an infinite 
sequence £n € Xs so that S”xn =< (i.e., x has infinitely many S-preimages). 


EXERCISE 5. Prove Proposition 4.9. 


COROLLARY 4.10. For Penrose tilings, the finite type tiling spaces are the same 
as the substitution tiling spaces: Xp = Xs and Xp, = Xs. 


Proor. By Proposition 4.9 and Corollary 4.8 it suffices to show S~! and oa 
exist on Xp and Xp, respectively. This is easy to see by inspection for Xp,. 
Although it is a little harder to see for Xp, it is also true in that case. 

Alternatively, this result for Xp also follows from Xp, using Example 5.20 
below. 
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5. Applications of topological dynamics 


5.1. Repetitivity and minimality. Tilings x constructed from aperiodic 
sets Q of prototiles, as well as substitution tilings, tend to have the following curious 


property. 


DEFINITION 5.1. A tiling x is called repetitive if for any patch y in x there is 
an r > 0 such that for any t € R? there is a translation TSy of y in x such that 
supp(T*y) C B, +t. 


In other words, a copy of y occurs “nearby” any given location t in x. Since 
all periodic tilings are repetitive, we think of repetitivity as a generalization of 
periodicity. 

Let (X, T) be a dynamical system. Let U C X be open and let x € X. Define 
the return set of x to U to be 


(5.1) R(z,U) = {t € R! : Ttz € U}. 


A set R C Rê is called relatively dense if there is an r > 0 such that every r-ball in 
R intersects R. 


DEFINITION 5.2. A point x € X is almost periodic if R(x, U) is relatively dense 
for every open U C X with R(x, U) #0. 


For a tiling space X C Xz, let y € T* and let R = supp(y). Define X (y) = 
{x € X : x|R] = y}. For e > 0 define the cylinder set 


Uy = TP”: X (y) = {T*x : x € X(y),t € Be}. 


Clearly Uye is open. 

Without loss of generality, we can assume by translating that the support of 
each patch y € T* contains the largest possible ball B, around the origin. It follows 
that the cylinder sets {Uy e„ :y E T*, €n — 0} form a basis for the tiling topology 
on X. In Definition 5.2, it suffices to check only the sets U belonging to this basis. 
Thus we have the following. 


PROPOSITION 5.3. Let (X,T) be a tiling dynamical system. Then a tiling x € X 


is repetitive if and only if it is an almost periodic point". 


A dynamical system (X,T) is called minimal if there are no proper closed T- 
invariant subsets of X. For a point x € X we define its orbit by O(x) = {Ttx: t € 
Ri}, and its orbit closure O(x) C X to be the closure of O(x) in X. 


GOTTSCHALK’S THEOREM. ([Got44]) A dynamical system (O(x), T) is mini- 
mal if and only if x is almost periodic. 
It follows from the minimality of the dynamical system (O(x),T) that O(y) = 


O(«) for all y € O(a). In this case it follows from Gottschalk’s Theorem that y is 
almost periodic too. 


DEFINITION 5.4. Two repetitive tilings x,y are said to be locally isomorphic if 
O(a) = O(y). 


"Because of this, the author previously used the term “almost periodic tiling” to mean a 
repetitive tiling. 
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Geometrically, two locally isomorphic tilings x and y have exactly the same 
patches. Dynamically, local isomorphism means z and y belong to the same minimal 
tiling dynamical system. 

A dynamical system (X,T) is transitive if there is a single orbit: O(a) = X 
for all x € X. This is a special case of minimality. A tiling x is periodic if and 
only if O(x) = O(a). In this case (T, O(x)) is transitive, and O(x) = R?/T is a 
d-dimensional torus (e.g., Xs = Z? /RÎ). 

We say a tiling is properly repetitive if it is repetitive but not periodic. An 
easy application of Zorn’s lemma shows that every dynamical system (X,T) has a 
minimal T-invariant subset Y C X. It follows that every tiling space contains a 
repetitive tiling (this argument appears in[RW92]). Of course in general this tiling 
may be periodic, but if we know X is an aperiodic tiling space, then it must contain 
a properly repetitive tiling. 


EXERCISE 6. Starting with the Penrose tiles P, construct a new example of an 
aperiodic prototile set P’ such that not every tiling x € Xp: is repetitive. 

A dynamical system (X, T) is called topologically transitive if there exists x € X 
such that O(a) = X. Clearly transitive implies minimal which implies topologically 
transitive. In each case the converse is false. Exercise 7 shows that (Xs,,T) is 
topologically transitive but not minimal. Later we will show that the Penrose tiling 
dynamical system (Xp,7T) is minimal but not transitive. A minimal dynamical 
system which is not transitive is called properly minimal. 


EXERCISE 7. Show that the tiling dynamical system (Xs,,T) is topologically 
transitive but not minimal. What can you say about (XR,,,T)? (See Examples 2.5 
and 2.6.) 


REMARK 5.5. When x periodic, O(a) is a torus, i.e., a connected manifold. 
It turns out that this is an exceptional situation. One can show that if (X,T) is 
aperiodic and topologically transitive, then for each patch y, X (y) is homeomorphic 
to a Cantor set. Since for € sufficiently small, Uy e = T«X(y) is homeomorphic to 
B. x X(y), every point x € X has a neighborhood homeomorphic to a product of 
Rt and a Cantor set. Such a space is called a lamination. In particular, a tiling 
space X is almost never connected. 


5.2. The repetitivity of substitution tilings. 


PROPOSITION 5.6. Let Xs be the substitution tiling space corresponding to an 
invertible primitive tiling substitution S. Then any Xs is an aperiodic tiling space 


PROOF. Suppose Tto = x for some to Æ 0. Since C is invertible, (4.2) implies 
(5.2) TECO "g =O "g 


for all n > 1. Choose n so large that T*int(L"D)  int(L"D) Æ 0 for all D €T. 
Since for some D € T, L” D € zx, this contradicts (5.2). 


REMARK 5.7. In Proposition 5.6, invertibility is also necessary [S0198]. 


The following generalizes a well known result for discrete substitution dynam- 
ical systems (see [Que87]). 
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THEOREM 5.8. Let Xs be a substitution tiling space corresponding to an prim- 
itive tiling substitution S. Then any x E€ Xs is repetitive. Moreover, any x,y E Xs 
are locally isomorphic. In particular, (Xg,T) is minimal. 


COROLLARY 5.9. The tiling space Xg corresponding to invertible primitive 
tiling substitution S consists of properly repetitive tilings. 


PROOF. This follows from Gottschalk’s Theorem, Theorem 5.6 and Theorem 5.8. 


PROOF OF THEOREM 5.8. Assume without loss of generality A > 0 (otherwise 
replace S with SF). Let x € X and let y be a patch in x. Fix Do € T and choose 
k so large that y is a patch in S*~!({Do}). Then y is a patch in S*({D}) for all 
DET. 

Since x € X it follows from Proposition 4.9 that there exists £k € X so that 
SFrp = x. Let x}, = L £k € Xger. Note that CFx, = zx. 

Let s be the largest diameter of LD € LT and let r = 2s. It follows from 
the triangle inequality that for any D’ ~ LË D, if t € D’ then D' C B, +t. 

Thus any r-ball B, + t in R? contains a tile D’ € z/,, and the patch C*({D’}) 
in x, which has support D’, contains a sub-patch that is a copy of y. 


5.3. Self-affine tilings. Let S = LC be a tiling substitution on 7. A tiling 
zo € Xz is called a self-affine tiling with expansion map L if Sxo = xp. When 
L is a similarity, xo is called a self-similar tiling. Self-affine tilings play the same 
role in the theory of substitution tiling dynamical systems that fixed points play 
in the theory of discrete substitutions. In particular, self-affine tilings provide a 
alternative definition for substitution tiling spaces. 


THEOREM 5.10. If S is a primitive tiling substitution then there exists k > 0 
and zo € Xs such that SF £o = zo. 


PROOF. First, we assume without loss of generality that A > 0, since otherwise 
we can replace S with S* so that AF > 0. Fixing Dı € T, there is a translation of 
L~+(D,) in C(D,). By taking additional powers of S, if necessary, we can assure 
that L~1(D,) C int(supp(D)). Similarly, there is a sequence Dp of tiles equivalent 
to Dı so that for all k > 0, L~'(Dx) € C(Dp-1) and L~'(D,) C int(supp(Dz_-1)) E 
int(D,). Since L is expanding, there exists a unique c € int(D}) satisfying 


loc) 
ce |) LD. 
k=1 
The point c is called a control point. 

Let Do = D, — c. This is a tile with a control point at the origin. Define the 
sequence x, = S*({Do}) and use this sequence to construct the substitution tiling 
space Xs. Because of the location of the control point, we have that x,_1 is a 
sub-patch of x, for all k, and £k = Saxpz_1. 

Let xo = Uk>1 Tk, and note that xo is a tiling of RI by the choice of an interior 
control point. By Proposition 4.9 we have zo € Xs, and also Szo = £o. 


It follows from Gottschalk’s Theorem and Theorem 5.8 that for any tiling sub- 
stitution S, O(a) = Xs for any x € Xs. By Theorem 5.10 it follows that there 
exists a self-affine tiling zo € Xs, (i.e., S*a9 = zo). Thus for any tiling substitu- 


tion S there exists a self-affine tiling zo such that Xs = O(xo). This hints at the 
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alternative definition of Xs mentioned above. However, as the next example shows 
it is not completely straightforward. 


EXAMPLE 5.11. Let S be the table substitution and let yı be the patch con- 
sisting of two rows of two horizontal table tiles, arranged in a 4 x 2 rectangle. Note 
that yı € F, the set of forbidden patches for table tilings. Put yn = $?—y,, so 
that yn C Yn+1. Then yo = Ux>i1yx satisfies S?yo = yo. However yo ¢ Xs since it 
contains the patch yı C yo and yı E F. 

One can show that the patch yı in the occurs only at the origin in yo, but 
nowhere else. In particular, the self-similar tiling yo is not repetitive. 


The following result shows the correct way to define a substitution tiling space 
in terms of self-affine tilings. 


PROPOSITION 5.12. Let S = LC be a primitive tiling substitution and let xo € 
Xr be repetitive and satisfy S*xo = xo for some k > 1. Then Xs = O(a). 


REMARK 5.13. Suppose we start with a tiling substitution S that we do not 
assume a priori satisfies the 2-patch closure property, but that satisfies S*a 9 = xo 
for some zo E€ Xz. It then follows that S does satisfy 2-patch closure. Once this 
is known, we can use Proposition 5.12 to define Xs provided we can verify that xo 
is repetitive. This may be easier in practice than verifying that S' satisfies 2-patch 
closure. 


Next, we consider the question of what linear maps L can occur as the expansion 
map for a self-affine tiling. 

Let S = AC be a primitive tiling substitution for d = 1, and let Sao = zo be 
a self-similar tiling of R. If A is the structure matrix for S, then it is easy to see 
that A must be the Perron-Frobenius eigenvalue of A: the unique real eigenvalue of 
largest modulus (see Section 6.3 below). 


EXAMPLE 5.14. A discrete substitution o is a mapping from a finite alphabet 
A = {1,...,n} to the set of all non-trivial finite words in the alphabet. For a 
concrete example, see (8.2) in Section 8.6 below (see also [Que87]). 

We define the structure matrix A of o to be the nxn matrix such that the entry 
Ai j is the number of times the letter i occurs in o(j). Assuming A is primitive, we 
let à > 0 and a = (a1,...,@n) > 0 be its Perron-Frobenius eigenvalue and eigen- 
vector. Given a primitive non-negative integer matrix A, it is easy to manufacture 
a substitution ø with structure matrix A (i.e., this amounts to choosing orders for 
letters in the words corresponding the columns of A). 

Now let 7 = {[0,a;] : i = 1,...,n} be a set of prototiles in R and for conve- 
nience, identify 7 with A by identifying [0,a;] with i. Define a tiling substitution 
S in such a way that S([0,a,;]) is the partition of \- [0,a,] into translates of the 
intervals corresponding to i1,12,...,im,, in order, where o(j) = i1i2...im,. It 
follows that S is a tiling substitution for d = 1 with expansion A. Thus we can 
construct the corresponding substitution tiling dynamical system (X 5,7), and by 
Theorem 5.10 there exists a self-similar tiling x9 € Xg with S*xo = xo for some k. 


Those numbers A that can be obtained as Perron-Frobenius eigenvalues of a 
primitive non-negative integer matrix where classified by Doug Lind [Lin84], who 
called them Perron numbers. They consist of all positive real algebraic integers 
such that any Galois conjugate A’ of À satisfies |\’| < A. Modulo some technicalities, 
we have essentially proved the following observation of Thurston. 
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PROPOSITION 5.15. (Thurston, [Thu89]) A positive real number A is the expan- 
sion for a self-similar tiling of R (or equivalently a 1-dimensional tiling substitution) 
if and only if it is a Perron number. 


A similar result holds for d = 2 in the the self-similar case. 


THEOREM 5.16. (Thurston [Thu89], Kenyon [Ken96]) Given A € C there is a 
primitive 2-dimensional self-similar tiling substitution S = AC and a self similar- 
tiling xo E€ Xs with expansion àA if and only if A is a complex Perron number: 
an algebraic integer À such that any Galois conjugate ’ of A, except possibly the 
complex conjugate À, satisfies |X'| < |A]. 


We call an expansion L € Gi(d,R) a Perron expansion if its eigenvalues A = 
{\i,...,Aa}, written with multiplicity, satisfy the condition that for every A € A 
with multiplicity k, if A” is a Galois conjugate of A with |A| > |A|, then à € A with 
multiplicity k’ > k. This idea generalizes both real and complex Perron numbers 
viewed as expansions of R and C & R? respectively. 


THEOREM 5.17. (Kenyon, [Ken90]) If a diagonalizable linear map L € Gl(d, R) 
is the expansion for a primitive self-similar tiling substitution S = LC, and Sap = 
xo for some xp E Xs, then L is a Perron expansion. 


Kenyon [Ken90] claims that the converse is also true. 


5.4. Local mappings. Local mappings play much the same role in tiling dy- 
namical systems that sliding block codes play in symbolic dynamics. The following 
version of the definition comes from [PSO1]. 


DEFINITION 5.18. A continuous mapping between tiling spaces Q : X — Y is 
called a local mapping if there is an r > 0 so that for all x € X, Q(x)[{0}] depends 
only on z[B,]. We say Q(x) is locally derivable from x. If Q is invertible, we say 
x and Q(x) are mutually locally derivable. 


EXERCISE 8. Show that a local mapping is continuous and T-equivariant (i.e., 
T'Qz = QT*tx). Moreover, if a local mapping is invertible then its inverse is also a 
local mapping. Thus a composition mapping is local. 


Now consider two dynamical systems (X, T) and (Y, T). A surjective continuous 
mapping Q : X — Y so that TtQ = QT* for all t € R? is called a factor mapping, 
and (Y,T) is called a factor of (X,T). An invertible factor mapping Q is called 
a topological conjugacy. In this case the two dynamical systems are said to be 
topologically conjugate. 


LEMMA 5.19. IfQ:X —Y is a surjective local mapping between tiling spaces 
then it is a factor mapping between tiling dynamical systems. If Q is invertible then 
it is a topological conjugacy. 


EXAMPLE 5.20. (Equivalence of different Penrose tilings). Let Xp denote 
rhombic Penrose tilings and let let Xp, be the triangular Penrose tilings. Fig- 
ure 12 shows how these two types of tilings are mutually locally derivable. Similar 
invertible local mappings connect these examples to several other famous types of 
Penrose tilings not discussed here (e.g., the “kites and darts” tilings and Penrose’s 
original “pentagon” tilings: see [GS87]). It follows that all the various correspond- 
ing Penrose tiling dynamical systems are topologically conjugate. 


TILINGS 19 


Ld a 


FIGURE 12. The local equivalence between rhombic Penrose tilings 
and triangular Penrose tilings. 


DEFINITION 5.21. A factor mapping Q : X — Y is almost 1:1 if there is a point 
yo E€ Y so that cardQ7!({yo}) = 1. 


The name of this kind of factor comes from the fact that if (Y, T) is topologically 
transitive, then {y € Y : card(Q7!(y)) = 1} is a dense Gs set. 

Now suppose T4 is a marked version of a prototile set T (e.g., in the way P is 
a marked version of Rs). Then the local mapping F that erases or “forgets” the 
markings is an example of a local mapping. 


THEOREM 5.22. (Goodman-Strauss [G-S98]) Suppose (Xs, T) is a self-similar 
substitution tiling dynamical system with Xs C Xz, for 2-dimensional set T of 
prototiles with finite local complexity. Then there exists a marking Ty of T, and a 
local matching rule Ta such that the forgetful mapping F : Xr, — Xr satisfies 
F(X7z,) = Xs, and F : Xz, — Xs is almost 1:1. 


In most cases this mapping is not invertible (i.e., is not a topological conjugacy). 
In the terminology of symbolic dynamics, a non finite type factor of a finite type 
tiling dynamical system should be called strictly sofic. It follows that a substitution 
tiling space is always sofic and usually strictly sofic. An exception is the following. 
Let F : P — Rs be the mapping that erases the arrows defining the Penrose 
matching rules. The tilings X = F(Xp) C Xr, are called the unmarked Penrose 
tilings. 


THEOREM 5.23. (de Bruijn [dB81]) The marked Penrose tilings are mutually 
locally derivable with the unmarked Penrose tilings. 


It is interesting to compare these examples to the the situation for 1-dimensional 
symbolic dynamics. In that case, the intersection between substitution dynamical 
systems and sofic shifts (or shifts of finite type) consists only of purely periodic 
examples. 


REMARK 5.24. Even though the finite type property is not closed under fac- 
torization, it is closed under topological conjugacy. See [RS01] for a proof. 


To illustrate one further application of local mappings, we return to the idea 
of an imperfect tiling substitution. 


DEFINITION 5.25. An imperfect tiling substitution is a mapping S = LC : 
Xr — Xz where L € Gi(d,Z) is expansive and C : Xr —> Xyz-17 is a local 
mapping. For emphasis, the corresponding tiling dynamical system Xs C Xz is 
called an imperfect substitution tiling dynamical system . 

Note that Examples 4.4 and 4.5 both satisfy this. An imperfect substitution 


is invertible if C is an invertible local mapping and it is also possible extend the 
notion of primitivity to this case. 
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FIGURE 13. Iterating the boundary in the binary tiling decompo- 
sition to obtain a perfect decomposition with a fractal boundary. 


THEOREM 5.26. (Priebe and Solomyak [PS01]) Suppose Xs C Xr is a imper- 
fect substitution tiling dynamical system with d = 2. Then there is new prototile 
set T’ and a perfect substitution S” on T” so that the substitution tiling dynamical 
system Xs: C Xz and Xs are topologically conjugate via a local mapping (i.e., 
corresponding tilings are mutually locally derivable). 


This theorem is proved using a construction called “iterating the boundary” 
which often results in producing fractal tiles. In Figure 13 this idea is illustrated 
in the case of the binary tiling system of Example 4.5 


REMARK 5.27. One way in which tiling dynamical systems differ from discrete 
symbolic dynamical systems is the following. The Curtis-Lyndon-Hedlund Theo- 
rem [Hed69] says that any factor mapping between discrete symbolic dynamical 
systems is implemented by a sliding block code. In tiling dynamical systems the 
equivalent question is whether every topological conjugacy is implemented by a 
local mapping. i.e., is the converse to Lemma 5.19 true? A negative answer was 
provided by Petersen [Pet99] and Radin and Sadun [RSO1]. 


5.5. Incongruent tilings. In the case that x is a periodic tiling one has 
O(a) = O(x). In other words there is a single orbit. It follows that, up to translation 
or congruence, there is just a single tiling. 

THEOREM 5.28. If x is a properly repetitive tiling then the number of orbits 
in O(a) is uncountable. There are uncountably many incongruent tilings in a local 
isomorphism class. 


This follows directly form the next lemma, which illustrates the power of simple 
topological ideas in this subject. 


PROPOSITION 5.29. Suppose (X,T) is a minimal dynamical system with T an 
action of R. Let Q C X be such that X can be expressed as a disjoint union of 
orbits X = UregO(«). Then either card(Q) = 1 or card(Q) > Xo. 


PROOF. If Q is finite then card(Q) = 1 since (X,T) is minimal. Thus we 
suppose Q is infinite. 

For x € Q, write x in terms of its tiles £ = {Dj, Do,...}. Let Vi(x) = {Ta : 
0 € TtD;}. Then O(#) = U%,Vi(x) is a countable decomposition of O(x) into 
nowhere dense sets (i.e., they have an empty interior). Thus X = Upeo UZ, Vi (2). 
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But since X is a compact metric space, it follows from the Baire Category Theorem 
(see [Mun75]) that X is not a countable union of nowhere dense sets. This implies 
Q is uncountable. 


COROLLARY 5.30. (Penrose) In any Penrose tiling space X there are uncount- 
ably many incongruent Penrose tilings. 


5.6. Quasicrystallography. Let M (d) denote the set of all rigid motions of 
R¢ (i.e., the set of congruence transformations). Let O(d) be the subgroup of M (d) 
fixing the origin. Denote the subgroup of translations in M (d) by Rê. 

Suppose x is a periodic tiling. The following ideas are basic to mathematical 
crystallography (see [Sen95]). The space group or symmetry group of x is defined 
Gs ={M € M(d): Mx = x}. The translation group Iy is a normal subgroup, and 
the quotient H} = G,/T., called the point group is isomorphic to a finite subgroup 
of O(n). The Crystallographic Restriction is the theorem which says that in any 
dimension d, there are only finitely many possibilities for H,. In particular, for 
d=2no M € H, can have order 5. 

Now we sketch the outlines of a theory of quasicrystallography (see [Rob96b] 
for more details). For a tiling x € X define Gz,.x = {M _ E€ M(d): Mz € X}. 


LEMMA 5.31. [Rob96b] If (X,T) is minimal then Gr, x = Gy, x for all x,y € 
X. 


EXERCISE 9. Prove lemma 5.31. 


In the minimal case, we write Gx. It follows that Gx is a closed subgroup of 
M(d) containing Rt as a normal subgroup. We call Gx the quasisymmetry group 
and we call the quotient Hx = Gx /R? the quasicrystallographic point group. The 
algebraic situation is simpler than in the case of symmetry groups. One always has 
that Gx is a semi-direct product of R and Hx. In particular, Hx is isomorphic 
to a closed subgroup of O(d). 


PROPOSITION 5.32. For the Penrose tiling space Xp, Hx, = Dio (the dihedral 
group of order 20). For the Pinwheel tiling space X, Hx = O(2). Thus Hx contains 
the circle T as a subgroup. 


In fact, one can construct examples X so that Hx contains any finite order 
rotation. It follows that there is no Crystallographic Restriction for quasicrystals. 


6. Applications of ergodic theory 


6.1. Measures. Let (X,T) be a dynamical system. In this section we discuss 
the set M(X) of Borel probability measures in X. 

Without going into a lot of details (see for example [Wal82]) we mention that a 
measure u E€ M(X) is a function that assigns a number 0 < u(E) < 1 to a Borel set 
E C X. One way to interpret u is as a “probability law” in which u(E) measures 
the probability that a randomly chosen point x € X belongs to E. Borel sets can 
be complicated, but include all open sets, closed sets, and most of the other sets 
that typically arise in practice. We will assume all sets mentioned are Borel sets. 

The integral of a function with respect to a measure p is denoted fy f(a) dy. 
If we think of f as a random variable on X, then the integral is its expectation. 

Of particular interest to us will be T-invariant measures. These are the mea- 
sures that satisfy u(T*E) = (E) or f(Ttx) du = fy f(x) dp for all t € R?. We 
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denote the set of all invariant measures by M(X,T). One can show that always 
M(X,T) #0. 

An important feature of measure theory is that a measure u is completely 
determined by its values on a collection of sets smaller than the collection of Borel 
sets. For a tiling dynamical system (X, T), a measure u € M(X,T) is determined 
by its values on cylinder sets. One can show that there is a function uo : T* — R 
such that 


(6.1) p(Uy,e) = poly) Vol(B1)ef, 


for e sufficiently small. This generalizes a similar and well known result that holds 
for discrete symbolic dynamical systems. 


6.2. Unique ergodicity. An invariant measure is said to be ergodic if T'E = 
E for all t € R? implies u(E) = 0 or (E) = 1 (note the similarity to the idea of 
minimality). A dynamical system is called uniquely ergodic if M(X,T) = {u}. In 
this case u is always an ergodic measure. 


THEOREM 6.1. If S is a primitive tiling substitution then the corresponding 
tiling dynamical system (Xs, T) is uniquely ergodic. 


A similar result is well known in the case of discrete substitution systems (see 
[Que87]). We will prove Theorem 6.1 below. Later we will also discuss a dif- 
ferent kind of uniquely ergodic tiling dynamical system, but first we discuss the 
consequences of unique ergodicity. 

Uniquely ergodic dynamical systems satisfy the following especially strong ver- 
sion of the Ergodic Theorem. 


THEOREM 6.2. If a dynamical system (X,T) is uniquely ergodic then for all 

complex valued continuous functions f on X 
1 

6.2 lim —— Tx) dt = du, 
(6.2) 100 Vol(By) Jp, f(T"2) fre p 
where the expression on the left, viewed as a function of x, converges uniformly to 
the integral the right (a constant). Conversely, if for all continuous f and for all 
x, the limit in (6.2) exists, then (X,T) is uniquely ergodic. 


Now let us assume (X, T) is a minimal uniquely ergodic tiling dynamical system. 
Let x € X be a tiling and let y € T* be a patch that occurs in x. We know that 
the occurrences of y are relatively dense, but suppose we want a more quantitative 
description of this repetitivity. For simplicity we assume, as above, that supp(y) 
contains a maximal ball around the origin. Let P(x,y) = {t € R? : Tty C x}. Note 
that this is a subset of R(x, Uy e). Recall that the characteristic function of a set 


U is 
oe 1 ifzeu, 
XUE Jo ifzguU. 


We have in particular fẹ xu(x)du = (U). It turns out that even though char- 
acteristic functions of cylinder sets are not continuous, Theorem 6.2 still holds for 
them. 
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COROLLARY 6.3. A tiling dynamical system (X,T) is uniquely ergodic if and 
only if for any x E€ X and any y € T* the following limit exists: 


(6.3) card(B, N P(x, y)). 


i 1 
lim =———— 
t— oo Vol( B+) 
If (X,T) is uniquely ergodic, then the value of limit (6.3) is yoly). 


This result explains the combinatorial and geometric meaning of the unique 
invariant measure in the uniquely ergodic case: it determines the frequency of all 
the different tiling patches y in all the tilings x € X. 

In fact, a similar result holds under the weaker assumption that u € M(X,T) is 
just ergodic. We say x is generic for u if (6.3) holds. That is, every patch has a well 
defined frequency. It follows from the Birkhoff Ergodic Theorem (see [Wal82]) that 
a.e. x € X is generic. If (X,T) is not uniquely ergodic, different tilings will be 
generic for different ergodic measures. This will result in different patch frequencies 
for different tilings. Moreover, by Corollary 6.3 there will always be some tilings 
x € X so that the limit (6.3) diverges. In these tilings, certain patches will not 
have well defined frequencies. 


6.3. Perron Frobenius theory. In this section we present the main part of 
the proof of Theorem 6.2. We use a well known argument (see [Que87]) based on 
the Perron-Frobenius Theorem to show that in a substitution tiling space, every 
prototile occurs in every tiling with a well defined density. Unfortunately this is not 
quite enough to prove Theorem 6.2. In the next section we show how to generalize 
the idea of a higher block code from symbolic dynamics to tiling spaces, and how 
to use this idea to finish the argument. 


THEOREM 6.4 (Perron-Frobenius Theorem (see [Rue69])). Let A > 0 bea 
real square matrix with A > 0 for some k > 1. Then there is a simple positive 
eigenvalue w > 0 with w > |w'| for all other eigenvalues w’. Let a and b be the 
eigenvectors corresponding to w for A and AT. Then a,b > 0 and for any v € R? 
(6.4) im, Ate a 

n= w” 

The eigenvalue w > 0 and eigenvector a > 0 are called the Perron-Frobenius 

eigenvalue and Perron-Frobenius eigenvector of A . 


COROLLARY 6.5. let A be the structure matrix for a primitive tiling substitution 
S = LC and letw be the Perron-Frobenius eigenvalue of A. Then w = det(L). 


This is because both w and det(L) measure how the substitution S expands 


volumes. 


COROLLARY 6.6. Let X be a primitive substitution tiling space. For D; € T 
let D? = supp(S” Di). Then for any x € X and Di, D; E T the limit 
í 1 m 
exists. 
PROOF. Since w” = det(L”), we have Vol(D?) = w”Vol(D;). Write d;i = 
Vol(D;). Let A?; be the i, jth entry of A". Note that A}; is the number of tiles 
equivalent to D; that occur in 5”({D;}), so 


(6.5) Az; = card( D? N P(x, {D;})) 
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Since A?; = e; - (A"e;), it follows from (6.4) 


1 1 


where a = (aj,...,d@a) and b = (61,..., ba). 


With a bit more care one can prove the following refinement of the above. 


COROLLARY 6.7. Let X be a primitive substitution tiling space. Then for any 
xEX andDeET the limit 


: 1 
(6.6) Jim voy N P(a, {D})) 


exists. 


6.4. Higher patch tiles. This section discusses a technical result needed to 
complete the proof of Theorem 6.1. 

Let (X,T) be a tiling dynamical system X C Xz. Fixr > Oand let y1, yo,---,Ym 
be the equivalence classes of patches z[B,], x € X. Let 


B= eR Tals aa, 


where we assume F; has nonempty interior. Define a new tiling x’ by subdividing 
each tile D € x into the smaller tiles DN E;, i = 1,...,m. Up to equivalence there 
are only finitely prototiles in x’ for all x € X and we denote this new prototile set by 
T,. Moreover, the mapping H, : Xr — Xz, is clearly an invertible decomposition 
mapping. It is called higher patch mapping. This is similar to the idea of a higher 
block code in symbolic dynamics (see [LM95], [Que87]). 


PROPOSITION 6.8. If S is a tiling substitution on Xr then Sy = H7+SH, is a 
tiling substitution on Xz. If S is primitive then so is Sy. 


Theorem 6.1 now follows by applying Proposition 6.8 and Corollary 6.7. 


7. Mixing properties 


7.1. Mixing and eigenvalues: the geometric interpretation. Let (X,T) 
be a dynamical system and u € M(X,T). A complex function f € L? (X, 2) is called 
an eigenfunction if there exists a corresponding eigenvalue w € R? such that 


(7.1) fie ere ee F(x), 


for u a.e. x. Note that this “eigenvalue” is actually a vector! In Physics this might 
be called a “wave vector.” 

A constant function f is always an eigenfunction, corresponding to w = 0. 
Ergodicity is equivalent to w = 0 being a simple eigenvalue. Moreover, in the 
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ergodic case all the eigenvalues are simple, and the set X of eigenvalues is a countable 
subgroup of R? (see [Wal82)}). 

If the only eigenfunctions are constants, then T is said to be weakly mixing. The 
opposite situation is called pure discrete spectrum; it occurs when the eigenfunctions 
have a dense span in L?(X,). An eigenfunction f is continuous if it is equal p a.e. 
to a continuous function. 

A dynamical system is called strongly mixing (or just mixing)if for any Borel 
sets A and B 
(7.2) Feiss u(T*AN B) = u(A)u(B). 

—0o 
A well known theorem (see [Wal82]) says that weak mixing is equivalent to (7.2) 
holding except on a set of t of density zero (this set depends on A and B). 

Now consider a tiling dynamical system (X, T). As we will discuss below, the 
eigenvalues are related to the the “diffraction” properties of tilings x € X. Heuristi- 
cally, such diffraction is caused by constructive reinforcement of waves reflecting off 
atoms, usually thought of as being located at the vertices of a tiling. When a tiling 
exhibits diffraction it can be interpreted as evidence that the tiling has some sort 
of long range spatial order in the arrangement of its tiles. Periodic tilings always 
diffract, but as we will see below, so do some properly repetitive tilings. 

Conversely, if a tiling system satisfies a mixing property (i.e., a lack of diffrac- 
tion) then it indicates that its tilings enjoy some sort of long-range spatial disorder. 
Consider, for example, (7.2) applied to a pair of cylinder sets Uey, and Uey in a 
mixing tiling dynamical system (X,7). For a randomly chosen x € X and for t 
sufficiently large, the probability of seeing yı and Tty2, (up to an e-translation), 
is approximately e° Vol( B1)? uo(y1)uo(y2). Thus, the knowledge that yı sits at one 
place in x is approximately statistically independent of the knowledge that a copy 
of yə sits at any particular distant location. 


7.2. Weakly mixing tiling spaces. There are two known mechanisms for 
producing weakly mixing tiling dynamical systems. The first is related to the 
algebraic properties of the eigenvalues of the expansion. It generalizes ideas from 
the theory of discrete 1-dimensional substitutions. 

Let DET and z € X7. Define 


E(x) = {t € R? : 3D, Dz € x, D2 = Dı — t}. 


If (X,T) is a properly minimal tiling dynamical system, then E(x) is the same for 
all tilings x € X, and we write E(X). In addition, =(X) satisfies 


(7.3) f¢/\lt|| : t € E(X), t 40} = S7! CR? 
(see [S0197], Proof of Theorem 4.4). Note that this is the case when X = Xs is a 
primitive substitution tiling space. 

THEOREM 7.1. (Solomyak [S0197]) A number w € R¢ is an eigenvalue for an 


invertible primitive self-affine substitution tiling system (Xs,T) with S = LC if 
and only if 


(7.4) lim e?™i(L"tw) — ] 


T00 


for allt € E(Xs). Moreover, the eigenfunctions can always be chosen to be contin- 
uous. 
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This is a combination of Theorems 4.3 and 5.1 in [Sol97]. It generalizes a 
similar result of Host [Hos86] for 1-dimensional discrete substitutions. 


REMARK 7.2. In certain non-weakly mixing cases Solomyak [S0197] describes 
some eigenvalues explicitly. 


Our interest here is on weak mixing, and for this purpose the beauty of Theo- 
rem 7.1 is that it reduces the question to number theory. A real algebraic integer 
à > 1 is called a Pisot number if all of its Galois conjugates A’ satisfy |\’| < 1. 
Pisot’s Theorem says that a positive real algebraic number À satisfying e27™>" 
for some for some w € R must be a Pisot number (see [Sal63]). A complex Pisot 
number is a complex algebraic integer À all of whose Galois conjugates A’, except 
its complex conjugate, satisfy |A| < 1. The following generalization of this idea is 
due to Mauduit [Mau89]. 


—1 


DEFINITION 7.3. A set A’ = {A1,..., Aq } of distinct algebraic integers with 
|A;| > 1 is called a Pisot family if X’ is Galois conjugate of some A € A’ with X ¢ A’, 
then |X’| < 1. Otherwise A’ is called non-Pisot. 


A real Pisot number by itself is a Pisot family, as is a complex Pisot number 
together with its complex conjugate. The next result generalizes Pisot’s Theorem. 


THEOREM 7.4. (Mauduit, [Mau89]) Jf A’ = {\i,...,Aa} is set of distinct 
algebraic numbers such that 


od’ n 
lim e Xiv =] 
n— oo 


for somé (v1,... va) € (C0), then A’ is a Pisot family. 


A nonempty set A’ of distinct algebraic integers can be written as a disjoint 
union 


(7.5) A’ =A, UAgU=+-U Ng 


where for each i there exists a monic irreducible polynomial p; € Z[t] (the minimal 
polynomial) such that p;(A) = 0 for all A € Ag. 


DEFINITION 7.5. We say A’ is totally non-Pisot if each A; in the decomposition 
(7.5) is non-Pisot. 


It is clear that a totally non-Pisot family A’ is non-Pisot, and moreover, any 
nonempty subset A” C A’ is totally non-Pisot. We call an expansion L totally non- 
Pisot if its set A’ of eigenvalues, written without multiplicity, is a totally non-Pisot 
family. 


THEOREM 7.6. Suppose S = LC is a primitive invertible tiling substitution such 
that the expansion L is diagonalizable and totally non-Pisot. Then the substitution 
tiling system (Xs,T) is weakly mixing. 


This is essentially due to Solomyak [S0197], although our formulation is differ- 
ent. 


8Mauduit [Mau89] proves this for (R\{0}), but the proof works in the complex case. 
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PROOF. Let A’ = {\j,..., Aw} be the eigenvalues of L, written without multi- 
plicity. Then A’ is totally non-Pisot. Given t € Rt express t = Sa P;t where Pj 
is the projection to the eigenspace for Aj, parallel to the all the other eigenspaces, 
and tj = Pjt € C’. Then L't = 0", X"ty. 

Now suppose w Æ 0 is an eigenvalue for (Xs,T). Then by Theorem 7.1 we 
have 
(7.6) 1= lim e2Ti(L”tw) — Jim ep aa AR bw) — lim eres A(t Pw) 
where (t, w) = t - W is the inner product on C%, and P; is the complex-conjugate 
transpose of P;. 

Let A”(w) = {Aj € A’: P*(w) # 0}, and note that A”(w) # Ø since L is 
diagonalizable. By (7.3) there exists t € E(X) such that (t, Přw) # 0 for all j 
such that A; € A” (w). Applying (7.6), it follows from Theorem 7.4 that A” (w) is 
a Pisot family. But since A”(w) C A’, this contradicts the fact that A’ is totally 
non-Pisot 


REMARK 7.7. 
(1) Suppose d > 1 and let S = AC be a self-similar tiling substitution with 
A E€ R. If À is not real Pisot then (Xs,T) is weakly mixing. 
(2) Suppose d = 2 and S = AC, where A € C\R. If A is not complex Pisot, 
then (Xs,T) is weakly mixing. 
(3) In both cases above, for d = 2, the converse is also true, [S0197]. 
COROLLARY 7.8. The binary tiling dynamical system (Xg,,T) is weakly miz- 
ing. 
PROOF. We use the fact that C?(T) C A~?T*, A € R, so that S? = AC?. The 
3 


1 2 
the Perron-Frobenius eigenvalue for A, which is not real Pisot. 


structure matrix for S? is A?, where A = . Thus, the expansion A for S$? is 


The reader should compare the disordered appearance of binary tilings to the 
more regular appearance of the Penrose tilings. As we will see below, Penrose tiling 
dynamical systems have pure discrete spectrum. 


REMARK 7.9. Let S = AC be a tiling substitution with d = 2 and A € R\C. We 
show here that in order to establish that (Xs, T) is weakly mixing it is not sufficient 
(i) that |A| is not real Pisot, nor (ii) that |A|? = w (where w is the Perron-Frobenius 
eigenvector for A) is not real Pisot. 

For (i) consider p(t) = tt + t? — 1. This has complex Pisot root À = iy7, 

= 14% but |A| = v7 is not real Pisot. For (ii) consider q(t) = t? — t? + 10t — 5. 
This has a complex Pisot root A, but w = |A|? is a root of r(t) = t? — 10t? + 5t — 25, 
and so is not real Pisot. 

Finally, we observe that A complex Pisot implies that A is complex Perron. 
Thus by Theorem 5.16, in each case above, there exists a tiling substitution S with 
expansion À. The fact that the corresponding substitution tiling system (Xs, 7) is 
not weakly mixing follows from Remark 7.7, part 3. My thanks to the referee for 
providing these two examples (see also [S0199]). 


The second known mechanism responsible for producing weakly mixing tiling 
dynamical systems involves quasisymmetry. 
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PROPOSITION 7.10. [Rob96b] The group of eigenvalues Xx of a tiling dynam- 
ical system is invariant under the action of the quasicrystallographic point group 
Ax. 


COROLLARY 7.11. (Radin, [Rad94]) The pinwheel tiling dynamical system is 
weakly mixing. 


PROOF. Let X be the pinwheel tiling space. Then T C Hx. Since Xx must 
be discrete and T-invariant, it follows that Ux = {0}. 


REMARK 7.12. Since the almost 1:1 extension in Theorem 5.22 is always a 
metric isomorphism, we can obtain examples (starting e.g. with the binary tilings) 
of minimal uniquely ergodic finite type tiling spaces that are weakly mixing. On 
the other hand, the next result shows that none of these examples can be strongly 
mixing. 


THEOREM 7.13. (Solomyak, [Sol97]). No self-affine substitution tiling dynam- 
ical system (X,T) can be strongly mixing. 


7.3. Diffraction. X-ray diffraction experiments provide a powerful method 
for studying the microscopic structure of solids. In particular, quasicrystals were 
discovered (see [SBGC84]) as a result of the observations of unusual diffraction 
patterns. 

Mathematically we model diffraction as follows. A Delone set is a uniformly 
discrete and relatively dense subset of R? (see [Sen95]). Starting with a tiling z, 
we let z = v(x) be the Delone set of all its vertex points. We place an “atom” ôt, 
consisting of a unit atomic measure, at each vertex point t € z. The corresponding 
“solid” consists of the measure 


(7.7) Hzr = Xo Ôt. 


tEv(x) 


The physical diffraction pattern corresponds mathematically to the locations X/, 
of the point masses in the Fourier transform siz. There are many technicalities 
involving how this is actually computed. For example, one should really compute 
the transform of the autocorrelation measure, which is supported on the difference 
set z — z of z. We will not discuss these issues here, but rather refer the reader to 
[Hof97] for a nice overview. The main result in this area is the following. 


THEOREM 7.14 ([Dwo93], [Hof97]). Let (X,T) be a uniquely ergodic tiling 
space with eigenvalue set 4. Then for anys € X, Xi, CÈ. 


Because of this result, physicists have been especially interested in tiling dynam- 
ical systems with pure discrete spectrum. On the other hand, one would expect 
a weakly mixing repetitive tiling (like the binary tiling) to have no spots in its 
diffraction pattern. In spite of this, it still seems reasonable to regard this example 
as a “quasicrystal”: as we will see in the next section, it has entropy zero. 


7.4. Entropy. 


DEFINITION 7.15. Let (X,T) be a tiling dynamical system. For n > 0 the 
complexity c(n) of (X,T) is the number of different equivalence classes of tilings 
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x|Bn] for all x € X. The topological entropy is defined as the exponential growth 
rate in complexity: 
(7.8) A(X) = lim + Vol(Bi)" In(c(n)). 
n> n 

Nonzero entropy is, in some sense, the ultimate indication of disorder in a dy- 
namical system. A tiling whose orbit closure is a positive-entropy tiling dynamical 
system should probably be considered too disordered to be regarded as a quasicrys- 
tal. Some of the full tiling shifts, discussed above, do have positive entropy. Since 
they are topologically transitive, there are tilings among them with a positive en- 
tropy orbit closure. However, as expected, most of the tiling dynamical systems 
that arise in the study of quasicrystals do indeed have zero entropy. 


THEOREM 7.16. Suppose X is either (a) a substitution tiling space for a prim- 
itive invertible tiling substitution, or (b) a finite type tiling space which is uniquely 
ergodic. Then h(X) = 0. 


Part (b) is due to Radin [Rad91] in the Z? case and generalized to the case of 
tiling dynamical systems by Shieh [Sh]. Part (a) follows from the next theorem. 


THEOREM 7.17. [HRO2] Let (Xs,T) be a substitution tiling space, where S = 
LC is primitive and invertible. Suppose L has eigenvalues \1,...,q where |ra| < 
|A;| for all i. Let 


_ log|det(Z)| _ log(|Ai|- +--+ |Aal) 
o logal) log |Aq| ` 
Then the complexity satisfies 
c(n) < K- n°. 
for some K > 0. In the self-similar case, |A| = -+-+ =|Aal, so c= d. 


The proof follows the dissertation of Clifford Hansen ([Han00]), who studied 
the case of discrete multi-dimensional substitutions. It is based on the following 
consequence of local finiteness. 


LEMMA 7.18. Let Xr be a finite local complexity tiling space. Given m > 1 
there exists a constant J = J(m) > 0 so that for all n sufficiently large 


(7.9) #{x € T™ : xC y for some y E€ T™} < Jen. 


PROOF OF THEOREM 7.17. For D; € T and v = (1,1,...,1)' we have #(C?D,;) = 
(A?v);, where A is the structure matrix for S. Since S is primitive, the Perron- 
Frobenius Theorem implies 


where w > 0 and a,b > 0. Since also v > 0, we have r > 0. Thus there exists N so 
that for all sufficiently large p, 


p oe 
(7.10) max #(C D)<N-w?. 


For p > 0, call a patch y € L?T* a p-basic patch if for some D € y, each 
D' € y satisfies D'N D Æ Ø. We denote the p-basic patches, up to equivalence, by 
Yi» -- -Ym where M is independent of p. Then M’ = max{#(y}) : j =1,..., M} 
is also independent of p. 
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For a p-basic patch yf, we have C?yf € T. Let 
(7.11) My = max{#(C?y?) : j =1,...,M}. 
Since #(C?y?) < M’- maxper #(C? D), it follows from (7.10) that 
(7.12) Mp < M'N o. 


Let 6 = max{diam(D) : D € T} and let € be the maximum of all r > 0 such 
that all D € T satisfy B, C D — s for some s € R°. 
Fix q > (+ 1)/e. Then for all n > 1, egn > n+ ô. Define 


_ loglan) 
log |Aal 


n 


It follows that 
(7.13) €|Aa\?” = eqn > n+ ô. 

By the choice of 6, Bn+5 2 supp(æ|Bn]) for all n > 0 and all z € Xs C Xz. 
Since S' is invertible, one can define a super-tiling CPx € Xzr7 for any p > 1 and 
x E€ Xs. It follows from (7.13) that LPB: D Bn+s for all p > pn. Thus for each 
x € Xs, the patch (C~?z)[Bn+5] is a sub-patch of some p-basic patch Tty? in the 
super-tiling CPx. Hence, 

[Bn] © 2[Bss] = C?(C7Pa[Bnss]) © T°C? 

Now we apply Lemma 7.18, (7.11) and (7.12) to conclude that 

e(n) < J: Mp S JM'N -w? del e. gP 


for all p > pn once n is sufficiently large. In particular, we take n large enough that 
(7.10) holds for p = pn. It follows that 


c(n) 


log(qn) 
K' wr = K'w Tog Aql 


IA 


7 log(w) log(qn) 
= K'e wea 
= K'¢€ . no 
= K-n° 
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where K = K'q® = JM'Nd. 


We conclude with two open problems. 
e Can a uniquely ergodic finite type tiling space X be strongly mixing? 
e Is there an example of a uniquely ergodic, finite type tiling space that not 
an almost 1:1 extension of a substitution tiling space (or, more generally, 
that is not a tiling space having some other type of hierarchical structure)? 


8. Quasiperiodic tilings and models 


In this section we describe the projection method, which is the third general 
method (after local matching rules and tiling substitutions) for constructing aperi- 
odic tilings. This method also provides an algebraic/geometric model for the tiling 
spaces it produces. After describing the projection method, we briefly describe 
some other geometric models. 
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8.1. Kronecker dynamical systems. Let G be a compact abelian group, 
written additively, with normalized Haar measure y. Suppose there exists an in- 
jective continuous homomorphism v : R? — G. The continuous y-preserving R4 
action 


(8.1) T'g=ut)+¢g 
on G is called a Kronecker dynamical system. Every point in such a dynamical 
system is almost periodic. Any coset W = (R1) +g is T-invariant. It follows that, 


up to topological conjugacy, (W,T) is a minimal Kronecker system. Any minimal 
Kronecker system is uniquely ergodic (see [Wal82]). 


REMARK 8.1. By an appropriate choice of s one can ensure that (W,T) is 
properly minimal. 


Every Kronecker system has pure discrete spectrum. The eigenfunctions are 
the characters x € G of G. In particular, for x € G we have y(i(t)) € Ri = R4, 
so there exists w € R? such that y(u(t)) =e"), and the set © of eigenvalues of 
(W,T) is the set of all such w (see [Wal82] for details). 

Two dynamical systems (X, 7) and (Y,T), preserving measures p and v respec- 
tively, are said to be metrically isomorphic (see [Wal82]) if there exist T-invariant 
subsets Xo C X and Yo C Y with u(Xo) = v(Yo) = 1, and an invertible measure 
preserving Borel mapping Q : Xo — Yo so that QT = TQ. Metric isomorphism is 
the primary notion of isomorphism studied in ergodic theory. If (X,T) and (Y, S) 
are uniquely ergodic then topological conjugacy implies metric isomorphism. Note 
that any two metrically isomorphic ergodic dynamical systems must have the same 
eigenvalues. 


HALMOS-vON NEUMANN THEOREM. (See [Wal82]) Any dynamical system 
with pure discrete spectrum is metrically isomorphic to a Kronecker system. Every 
countable subgroup © C R? is the eigenvalue group for a Kronecker system with RÎ 
acting. 


REMARK 8.2. Every Kronecker system has entropy zero (see [Wal82]). 


8.2. The projection method. A landmark in the theory of aperiodic tilings 
is de Bruijn’s algebraic theory of Penrose tilings [dB81]. Originally, this theory 
described Penrose tilings as being dual (in the sense of graph theory) to so-called 
“grid” tilings. The generalization of this idea is called the grid method, and the 
tilings it produces are called quasiperiodic tilings. There are two alternate equivalent 
constructions of quasiperiodic tilings, mostly developed by physicists: the projection 
method and the cut method (see [ODK88]). Here we discuss the projection method 
because it is conceptually the simplest. 

Let Ell be a d-dimensional subspace of R” and let ų : R” — Ell be an isometric 
isomorphism. Let E+ be the perpendicular subspace, so that R” = Ell @ Et. 
Denote the projections to these two subspaces by m! and m+. Consider the integer 
lattice Z” C R”. Let Wo be the closure of (R?) in T” = R"/Z” and let Wo + g, 
g € EŁ, be an arbitrary coset. The tiling systems we will construct is closely 
related to the Kronecker system (Wo, T). 

Let K C E+ be compact with a Lebesgue measure zero boundary and a 
nonempty interior. Let Sx = K + Ell. For s € R let 


zs = (trl) (Sp N (Z” +8)) C RÊ. 
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Note that zs = zs if s—s’ € Z. Thus we can index using s € T” = R” /Z”, and 
this makes the mapping s > zs 1:1. 

The set zs is a Delone set. It is possible to topologize the collection of all 
Delone sets (using something similar to a tiling metric) in such a way that the set 
{zs : s € T”} is homeomorphic to T”. We also have zrts = T*zg, so in fact s = Zs 
is a topological conjugacy. 

We call s € T” regular if O(SK)N(Z" +s) =. Let Tj denote the set of all 
regular points in T”, and note that Tg has full Lebesgue measure. In many cases 
it will turn out that the points in zs are the vertices of a tiling x of R?. Let us now 
specialize to such a situation. 

Let d = 2 and consider the full tiling space Xr„, n > 4, from Example 2.6. 
Let K = r+ (Q) where Q = {q € R” : 0 < q; < 1} is the unit cube. Let s be as in 
Example 2.6, and let B be matrix having the vectors vz, as rows. Define s(t) = Bt. 


PROPOSITION 8.3 (de Bruijn, [dB81]). For each s € TẸ there exists a tiling 
xz E€ Xr, with v(x) = zs, i.e., x has the Delone set zs as its vertex points. 


In addition, de Bruijn [dB81] showed that the non-regular points in T” cor- 
respond to more than one tiling (but always a finite number). Moreover, each of 
these tilings can be obtained as a limit of the regular cases. This result can be 
made into a statement about tiling dynamical systems as follows: 


THEOREM 8.4. [Rob96b] Let W = (R?) +g forge Wt. Fors € W AT} 
define x = H (s) € Xr, to be the tiling with vertex set u(x) = zs. Then 


X =H(WNT?) 


is a tiling space that is minimal and uniquely ergodic. Moreover, H! extends to 
a continuous mapping P : X — W which is an almost 1:1 factor mapping and a 
metric isomorphism. 


In the case where + is chosen so that (W,T) is properly minimal, each tiling 
x € X will be aperiodic (i.e., properly repetitive). The tilings x € X are called 
quasiperiodic tilings and (X,T) is called a quasiperiodic tiling dynamical system. 


COROLLARY 8.5. Every quasiperiodic tiling dynamical system (X,T) is prop- 
erly minimal, uniquely ergodic and contains uncountably many incongruent tilings. 
Moreover, it has an almost 1:1 Kronecker system factor (W,T), to which it is met- 
rically isomorphic. It thus has pure discrete spectrum and entropy zero. 


It is possible to find explicitly the set X for these quasiperiodic tiling dynamical 
systems (see [Rob96a]). In particular, 


Y= BZ" = © NjVj inj E N}. 
j=l 


When d = 2 and n = 2” we have (R2) = T”. The case n = 4 gives a well 
known example called the octagonal or Ammann-Beenker tilings. Like the Penrose 
tilings, the Amman-Beenker tilings can also be generated by a local matching rule 
and by a tiling substitution (see [Sen95] for details). 

In the case that d = 2 and n = p is an odd prime, (R?) = TPt. Here E+ is 
the 1-dimensional subspace generated by (1,1,...,1), and T” = W œ W+ where 
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FIGURE 14. A patch of Ammann-Beenker tiling. 


Wt = E+/Z? = T. let us define y : T” — Wt by (s) = s1 + 52 +--+ + Sn 
mod 1. For 0 < t < 1, we define 


Xni = H(W +E) NTA) C Xr.. 


We can now state de Bruijn’s remarkable algebraic structure theorem for Penrose 
tilings. 


THEOREM 8.6. (de Bruijn [dB81]) X5, is the precisely the set F(Xp) of un- 
marked Penrose tilings. 


COROLLARY 8.7. [Rob96a] The Penrose tiling dynamical system (Xp,T) has 
pure discrete spectrum, with an almost 1:1 Kronecker factor T on the 4-torus T4. 
The eigenvalue group X, is the subgroup of R? is generated by the “Sth roots of 
unity” Vo,..., V4- 


REMARK 8.8. In a slight abuse of notation we write Xp = Z{e?™’/*] C R?. 


The tilings on which the factor map P : Xo, > T° fails to be 1:1 are a well 
known special class Penrose tilings. In particular, P is 2:1 on the Penrose tilings 
that contain infinite worms, and P is 10:1 on the various cartwheel tilings. In 
particular, P(x) = 0 for the cartwheel tilings x centered at the origin. See [GS87] 
for the geometric description of these special Penrose tilings. Similar results hold 
for all other classes of quasiperiodic tilings. 

The cases X54 for t Æ 0, are called generalized Penrose tilings. All the gen- 
eralized Penrose tiling dynamical systems have pure discrete spectrum and they 
all have the eigenvalue group © = Z{e?*’/>], which they share with the true Pen- 
rose tilings. Thus, by the Halmos von Neumann Theorem, all the corresponding 
tiling dynamical systems are metrically isomorphic. However, one can also show 
(see [Rob96a], [Le97]) that in some cases they are not topologically conjugate. 
This illustrates the interesting fact that ergodic theory and topological dynamics 
sometimes provide different invariants when applied to tiling theory. 
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EXERCISE 10. Show that the point groups satisfy Hx,, = Dio only if t = 0 or 
t = 1/2, and otherwise Hx, , = Ds. 


REMARK 8.9. A Sturmian shift dynamical system is a symbolic dynamical sys- 
tem obtained by coding an irrational rotation on the circle using a partition into two 
intervals (see [Ber00]). Consider the partition 7 of a quasiperiodic tiling space X 
according to which prototile contains the origin. The factor mapping P : X — W 
pushes 7 forward into a partition on W. We can recover the tilings by seeing how 
this partition tiles the orbits of R? in W. In this way, quasiperiodic tiling dynamical 
systems generalize the idea of Sturmian systems. 


8.3. Quasiperiodicity and the finite type property. An interesting gen- 
eral question concerns the relation between quasiperiodic tilings and the finite type 
property. Such questions were studied primarily by researchers interested in qua- 
sicrystals. See [Le97] for a good survey with complete references. Here, we discuss 
mainly the case of generalized Penrose tilings. 


THEOREM 8.10. [Le95] Let r = 444. 
(1) The tiling space X54 is a finite type tiling space if and only if t € Z[r]. 
(2) There exists a finite type tiling space Xz, by a marked version Ty of the 
tiles Rs so that the forgetful mapping F : Xt, — X54 is an almost 1:1 
factor if and only ift € Q(T). Ift € Q(r)\Z[r] then F is not a topological 
conjugacy; such examples are strictly sofic. 


REMARK 8.11. One should compare this with Theorem 5.22. 


REMARK 8.12. In general, a necessary condition for a quasiperiodic tiling sys- 
tem to be finite type or sofic is that the “slope” of Ell must be algebraic (in a certain 
precise sense). Thus sofic examples are very special (see [Le97]). 


REMARK 8.13. It seems to be still unknown exactly which quasiperiodic tiling 
systems are substitution tiling systems, and vice versa. 


8.4. Algebraic and geometric models. The subspace E+ in the definition 
of quasiperiodic tiling dynamical systems can be replaced by an arbitrary locally 
compact abelian group. In this case one obtains more general Delone sets called 
model sets. Such sets can be studied without the benefit of tilings or even dynamical 
systems theory. However, one can show that quite generally, the corresponding 
tiling dynamical systems have almost 1:1 Kronecker factors, defined on more general 
groups than tori. The chair tiling dynamical system can be obtained this way. (See 
[BMS98] for a discussion of this point of view.) 

Alternatively, one can ask, given a tiling dynamical system (X, T), whether 
there is a “classical” dynamical system that is its almost 1:1 factor. We refer to 
such a factor as an algebraic/geometric model for the tiling system. For example 
the Penrose tilings are modeled by T+. Here is another example. 

Let G be the locally compact group whose dual G is the group Z[5] @ Z[5] of 
dyadic rational vectors in R?. By Halmos-von Neumann theory, there is a unique 
minimal Kronecker R? action T on G with © = Z[$] @ Z[4] (this Kronecker system 
is called the R?-adding machine in [Rob99]). 


THEOREM 8.14. [Rob99] The chair tiling dynamical system (X s., T) has (G, T) 
as an almost 1:1 factor. Thus it has pure discrete spectrum with & = G. The table 
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tiling dynamical system (Xs,,T) has (G,T) as an almost 4:1 factor. The table 
system has precisely the same eigenvalues X as the chair, and thus is not weakly 
mixing. However, the table tiling system also does not have pure discrete spectrum; 
it has mixed spectrum. 


REMARK 8.15. The spectral properties of the table and chair were first com- 
puted in [S0197], where it was observed that the table system is similar to the well 
known discrete Morse sequence substitution dynamical system (see [Que87]). 


In the chair (or the table), the points where the factor mapping P fails either 
to be 1:1 (or 4:1) can be completely classified (see [Rob99]). In particular, chair 
tilings can have worms and cartwheels very similar to the ones that occur in Penrose 
tilings. The proofs in [Rob99] are completely general and show that similar results 
hold for all polyomino substitutions. In particular, there is an easy criterion for 
pure discrete spectrum. This turns out to be equivalent to the idea of “coincidence” 
in the theory of substitutions, and to the idea of “synchronizing” or “magic” words 
in the theory of discrete 1-dimensional finite type shifts (see [Rob99]). 


REMARK 8.16. Solomyak [S0197] defines a more general notion of coincidence, 
and uses it to prove that some examples (e.g., chair tiling systems) have pure 
discrete spectrum. 


Models based on Kronecker systems won’t work for weakly mixing examples 
because weak mixing implies the absence of nontrivial eigenvalues. Thus we ask, 
what could a model for a weakly mixing tiling dynamical system possibly look like? 


8.5. The dynamics of the substitution map and Markov partitions. 
Before describing an example of a model for a weakly mixing tiling dynamical sys- 
tem, we consider the dynamical properties of the action of an invertible substitution 
mapping S acting on a tiling space Xs. 


THEOREM 8.17. [Rob96a] The Penrose substitution S on the set Xp of Pen- 
rose tilings has the hyperbolic toral automorphism 


01 -1 0 
1 0 0 -1 
A=|ð 1 1 1 
1 0 -1 -1 


acting on T* as is an almost 1:1 factor. 


Similarly, one can show that the chair tiling dynamical system has a “hyper- 
bolic” automorphism of G as an almost 1:1 factor. These examples illustrate that 
tiling substitutions S tend to be hyperbolic. 

In [AP98] it is shown that the action S of a tiling substitution on Xg (i.e., the 
dynamical system (Xs, S)) is always a kind of generalized hyperbolic system called 
a Smale space. If the substitution S is based on a perfect decomposition, then 
then the tilings induce a partition on Xs (according to which prototile occurs at 
the origin). The perfect decomposition property implies this partition is a Markov 
partition (see [Bow78]). Conversely, if a Smale space has a Markov partition whose 
partition elements are connected, then the partition induces tilings on the stable 
manifolds, and these tilings are self-affine. 

Connectedness is needed because our definition of tilings requires connected 
tiles. Unfortunately, it is unknown whether every hyperbolic toral automorphism 
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has a Markov partition with connected partition elements. However, some partial 
results appear in [F198]. 


REMARK 8.18. It is known that for all Markov partitions for hyperbolic toral 
automorphisms of T? (see [Bow78]) and for typical hyperbolic toral automorphisms 
of T”, n > 3, (see [Caw91]) the boundaries of partition elements must be fractal. 
This implies that self-affine quasiperiodic tilings satisfying a perfect decomposition 
will almost always have tiles with fractal boundaries. 


8.6. A geometric model for a weakly mixing system. The example dis- 
cussed in this section is based on a kind of hyperbolic dynamical system J called 
a pseudo-Anosov diffeomorphism (see [FS79]). In this example, J is defined on a 
surface M of genus 2. Pseudo-Anosov diffeomorphisms always have Markov parti- 
tions, and we obtain self-similar tilings, as in the previous section, by intersecting 
the Markov elements partition with the unstable manifolds for J. Since for almost 
every point m € M, the stable manifold through m is homeomorphic to R, this 
example consists of 1-dimensional tilings. See [Fit98] or [FHR00] for details. 

Consider the discrete substitution ø given by 


1 > 1494 

2 — 142424 
(8.2) 3 — 14334 

4 — 1434. 


The structure matrix for this substitution is 


Re 


Nr OF 


which has the non-Pisot Perron-Frobenius eigenvalue À = +(74+V5+V2 V19+ 7/5). 
Following example 5.14, we construct the corresponding 1-dimensional tiling sub- 
stitution S with expansion A. It follows from Theorem 7.6 that (Xs,T) is weakly 
mixing. 


THEOREM 8.19. [FHRO00] Let Xs be the tiling space corresponding to the tiling 
substitution S in described above. Then (Xg,S) has an almost 1:1 factor that is a 
pseudo-Anosov diffeomorphism (M, J) on a surface M of genus 2. 


The corresponding tiling dynamical system (X 5,7) is metrically isomorphic 
a unit speed flow along the unstable manifolds of (M,J). Since not all the sta- 
ble manifolds are homeomorphic to R, this flow is not defined everywhere on M. 
However, after removing a “singular set” of measure zero, one can show that this 
flow is metrically isomorphic to a suspension of a self-inducing interval exchange 
transformation of four intervals (see [FHR00]). 


REMARK 8.20. One can obtain a 2-dimensional weakly mixing example by 
taking the Cartesian square. The geometric model for such an example is a 4- 
manifold. By Theorem 5.22, this 2-dimensional tiling dynamical system has the 
property that it is an almost 1:1 factor of a weakly mixing finite type tiling system. 


AP98] 
Brg66] 


Ber00] 


BMS98] 


Bow78] 
Caw91] 
Culg6] 
dB81] 
Dwo93] 
F879] 
Fit 98] 


FHR00] 


F198] 


GBN89] 


GL92] 
G-S98] 
Got44] 
GS87] 

Han00 
HRO2 
Hed69 
Hof97 
Hos86 
Kar96 
Ken90 


Ken92 
Ken93 


TILINGS 37 


References 


J. Anderson and I. Putnam, Topological invariants for substitution tilings and their 
associated C* algebras, Ergodic Thy. and Dynam. Sys. 18 (1998), no. 3, 509-537. 

R. Berger, Undecidability of the domino problem, Memoirs Amer. Math. Soc. 66 (1966), 
72pp. 

Valérie Berthé, Sequences of low complexity: automatic and Sturmian sequences, Top- 
ics in symbolic dynamics and applications (Temuco, 1997), Cambridge Univ. Press, 
Cambridge, 2000, pp. 1-34. 

Michael Baake, Robert V. Moody, and Martin Schlottmann, Limit-(quasi)periodic point 
sets as quasicrystals with p-adic internal spaces, J. Phys. A 31 (1998), no. 27, 5755- 
5765. 

Rufus Bowen, Markov partitions are not smooth, Proc. Amer. Math. Soc. 71 (1978), 
no. 1, 130-132. 

Elise Cawley, Smooth Markov partitions and toral automorphisms, Ergodic Theory Dy- 
nam. Systems 11 (1991), no. 4, 633-651. 

Karel Culik, II, An aperiodic set of 13 Wang tiles, Discrete Math. 160 (1996), no. 1-3, 
245-251. 

N. G. de Bruijn, Algebraic theory of Penrose’s nonperiodic tilings of the plane. I, I, 
Nederl. Akad. Wetensch. Indag. Math. 43 (1981), no. 1, 39-52, 53-66. 

Steven Dworkin, Spectral theory and X-ray diffraction, J. Math. Phys. 34 (1993), no. 7, 
2965-2967. 

A. Fathi and M. Shub, Some dynamics of pseudo-anosov diffeomorphisms, Traveaux de 
Thurston de les surfaces, Astersique, vol. 66-67, 1979, pp. 181-208. 

Thomas Fitzkee, Weakly mixing tiling flows arising from interval exchange transfor- 
mations, Ph.D. thesis, The George Washington University, 1998. 

Thomas Fitzkee, Kevin Hockett, and E. Arthur Robinson, Jr., A weakly mixing tiling 
dynamical system with a smooth model, Theoretical Computer Science, 303 (2003), 
447-462. 

Maki Furukado and Shunji Ito, The quasi-periodic tiling of the plane and Markov sub- 
shifts, Japan. J. Math. (N.S.) 24 (1998), no. 1, 1-42. 

F. Girault-Beauquier and M. Nivat, Tiling the plane with one tile, Topology and Cate- 
gory Theory in Computer Science, Oxford University Press, New York, 1989, pp. 291- 
333. 

C. Godréche and F. Lancon, A simple example of a non-pisot tiling with five-fold sym- 
metry, Journal de Physique 2 (1992), 207-220. 

Chaim Goodman-Strauss, Matching rules and substitution tilings, Ann. of Math. (2) 
147 (1998), no. 1, 181-223. 

W. H. Gottschalk, Orbit-closure decompositions and almost periodic properties, Bull. 
Amer. Math. Soc. 50 (1944), 915-919. 

Branko Griinbaum and G. C. Shephard, Tilings and patterns, W. H. Freeman and 
Company, New York, 1987. 

Clifford W. Hansen, Dynamics of multi-dimensional substitutions, Ph.D. thesis, The 
George Washington University, 2000. 

Clifford W. Hansen and E. Arthur Robinson Jr., On the complexity of self-affine tilings, 
preprint, 2003. 

G. A. Hedlund, Endormorphisms and automorphisms of the shift dynamical system, 
Math. Systems Theory 3 (1969), 320-375. 

A. Hof, Diffraction by aperiodic structures, The mathematics of long-range aperiodic 
order (Waterloo, ON, 1995), Kluwer Acad. Publ., Dordrecht, 1997, pp. 239-268. 

B. Host, Valeurs propres des systémes dynamiques définis par des substitutions de 
longueur variable, Ergodic Theory Dynam. Systems 6 (1986), no. 4, 529-540. 

Jarkko Kari, A small aperiodic set of Wang tiles, Discrete Math. 160 (1996), no. 1-3, 
259-264. 

Richard Kenyon, Self-similar tilings, Ph.D. thesis, Princeton University, 1990. 

Richard Kenyon, Rigidity of planar tilings, Invent. Math. 107 (1992), no. 3, 637-651. 
Richard Kenyon, Erratum: “Rigidity of planar tilings” [Invent. Math. 107 (1992), no. 
3, 687-651], Invent. Math. 112 (1993), no. 1, 223. 


38 


Ken96] 
Le95] 


Le97]| 


Lin84] 
LM95] 
Mau839] 
Mun75] 
ODK83] 
Pen74] 
Pet99] 
PSO1] 
Que87] 
Rad91] 
Rad94] 


Rad99] 
RW92] 


RS01] 
rRob71] 
Rob96a] 


Rob96b] 


Rob99] 
Rud88] 


Rud89)] 


Rue69)] 


Sal63] 


Sen95] 


E. ARTHUR ROBINSON, JR. 


Richard Kenyon, The construction of self-similar tilings, Geom. Funct. Anal. 6 (1996), 
no. 3, 471-488. 

Thang T. Q. Le, Local rules for pentagonal quasicrystals, Discrete Comput. Geom. 14 
(1995), 13-70. 

Thang T. Q. Le, Local rules for quasiperiodic tilings, The mathematics of long-range 
aperiodic order (Waterloo, ON, 1995), Kluwer Acad. Publ., Dordrecht, 1997, pp. 331- 
366. 

Douglas Lind, The entropies of topological Markov shifts and a related class of algebraic 
integers, Ergodic Theory Dynam. Systems 4 (1984), no. 2, 283-300. 

Douglas Lind and Brian Marcus, An introduction to symbolic dynamics and coding, 
Cambridge University Press, Cambridge, 1995. 

Christian Mauduit, Caractérisation des ensembles normauz substitutifs, Invent. Math. 
95 (1989), no. 1, 133-147. 

James R. Munkres, Topology: a first course, Prentice-Hall Inc., Englewood Cliffs, N.J., 
1975. 

C. Oguey, M. Duneau and A. Katz, A geometric approach to quasiperiodic tilings, 
Commun. Math. Phys. 118 (1988), 99-118. 

Roger Penrose, Pentaplexy, Bulletin of the Institute of Mathematics and its Applications 
10 (1974), 266-271. 

Karl Petersen, Factor maps between tiling dynamical systems, Forum Math. 11 (1999), 
no. 4, 503-512. 

N. Priebe and B. Solomyak, Characterization of planar pseudo-self-similar tilings, Dis- 
crete Comput. Geom. 26 (2001), no. 3, 289-306. 

Martine Queffélec, Substitution dynamical systems—spectral analysis, Lecture Notes in 
Mathematics, vol. 1294, Springer-Verlag, Berlin, 1987. 

Charles Radin, Disordered ground states of classical lattice models, Rev. Math. Phys. 
3 (1991), no. 2, 125-135. 

Charles Radin, The pinwheel tilings of the plane, Ann. of Math. (2) 139 (1994), no. 3, 
661-702. 

Charles Radin, Miles of tiles, American Mathematical Society, Providence, RI, 1999. 
Charles Radin and Mayhew Wolff, Space tilings and local isomorphism, Geom. Dedicata 
42 (1992), no. 3, 355-360. 

Charles Radin and Lorenzo Sadun, Isomorphism of hierarchical structures, Ergodic 
Theory Dynam. Systems 21 (2001), no. 4, 1239-1248. 

Raphael M. Robinson, Undecidability and nonperiodicity of tilings of the plane, Inven- 
tiones Math. 12 (1971), 177-209. 

E. Arthur Robinson, Jr., The dynamical properties of Penrose tilings, Trans. Amer. 
Math. Soc. 348 (1996), no. 11, 4447—4464. 

E. Arthur Robinson, Jr., The dynamical theory of tilings and quasicrystallography, Er- 
godic theory of Z? actions (Warwick, 1993-1994), London Math. Soc. Lecture Note Ser., 
vol. 228, Cambridge Univ. Press, Cambridge, 1996, pp. 451—473. 

E. Arthur Robinson, Jr., On the table and the chair, Indag. Math. (N.S.) 10 (1999), 
no. 4, 581-599. 

Daniel J. Rudolph, Rectangular tilings of R” and free R"-actions, Dynamical systems 
(College Park, MD, 1986-87), Springer, Berlin, 1988, pp. 653-688. 

Daniel J. Rudolph, Markov tilings of R” and representations of R” actions, Measure 
and measurable dynamics (Rochester, NY, 1987), Amer. Math. Soc., Providence, RI, 
1989, pp. 271-290. 

David Ruelle, Statistical mechanics: Rigorous results, W. A. Benjamin, Inc., New York- 
Amsterdam, 1969. 

Raphaél Salem, Algebraic numbers and Fourier analysis, D. C. Heath and Co., Boston, 
Mass., 1963. 


SBGC84] D. Schechtman, I. Blech, D. Gratias, and J. W. Cahn, Metallic phase with long range 


orientational order and no translational symmetry, Phys. Rev. Letters 53 (1984), 1951- 
1953. 

Marjorie Senechal, Quasicrystals and Geometry, Cambridge University Press, Cam- 
bridge, 1995. 


Sh] 
Sol97] 
Sol98] 


S0199] 


Thu89] 


Wal82] 
Wan61] 


TILINGS 39 


Jiunn-I Shieh, The entropy of uniquely ergodic tiling systems, preprint, University of 
Texas. 

Boris Solomyak, Dynamics of self-similar tilings, Ergodic Theory Dynam. Systems 17 
(1997), no. 3, 695-738. 

B. Solomyak, Nonperiodicity implies unique composition for self-similar translationally 
finite tilings, Discrete Comput. Geom. 20 (1998), no. 2, 265-279. 

Boris Solomyak, Corrections to: “Dynamics of self-similar tilings” [Ergodic Theory 
Dynam. Systems 17 (1997), no. 3, 695-738], Ergodic Theory Dynam. Systems 19 
(1999), no. 6, 1685. 

W. Thurston, Groups, tilings and finite state automata, AMS Colloquium Lecture 
Notes, Boulder, Colorado, 1989. 

Peter Walters, An introduction to ergodic theory, Springer-Verlag, New York, 1982. 

H. Wang, Proving theorems by pattern recognition. ii., Bell System Technical Journal 
40 (1961), 1-42. 


(E.A. ROBINSON) DEPARTMENT OF MATHEMATICS, GEORGE WASHINGTON UNIVERSITY, WASH- 
INGTON, DC 20052., USA 
E-mail address: robinson@gwu.edu 


